﻿L AMMERAL PRINCIPII ÎN MAȘINĂ GRAFICĂ PRINCIPII DE PROGRAMARE ÎN GRAFICA MAȘINILOR Principii de programare în grafica computerizată L Ammeral Hogeschool Utrecht Olanda JOHN WILEY & SONS Chichester • New York • Brisbane • Toronto • Singapore Seria COMPUTER GRAPHICS ÎN LIMBAJ C L Ammeral PRINCIPII PROGRAMARE ÎN GRAFICA MAȘINILOR Traducere din engleză de V A Lviv Moscova Sistemul solar BBC A UDC L Ammeral A Principii de programare în grafica computerizată Pe din engleza - M : "Sistem Sol", - p : ill ISBN - - -X (rusă) O introducere practică în grafica computerizată Sunt luate în considerare problemele de geometrie și programare analitică și proiectivă în grafica computerizată, problema formării unei B-spline, proiecțiile în perspectivă și eliminarea liniilor ascunse Algoritmii au fost aduși în programe de grafică "gata de lucru" în limbajul C Anexa oferă o scurtă prezentare a elementelor limbajului C Pentru o gamă largă de cititori care utilizează computere personale IBM PC sau compatibile cu acestea pentru a lucra cu informații grafice - > A I - u zhskh + au Pentru a obține valorile a, b, c, d, luați în considerare mai întâi punctul (x, y) = ( , ) Presupunând a= uy= în ecuația ( ), obținem x' = a y = C Dar în acest caz simplu, așa cum se poate vedea din fig (a), valorile lui x' și y' sunt cos cos^> sau folosind un vector coloană ( , ) 'x'I = Tcos^ -sin y'J L sin y? cos^> În cărțile de grafică pe computer, intrarea cu vector rând ( ) este mai frecventă decât cu vectorul coloană ( ) Aici vom folosi și o notație ca ( ) Într-o astfel de notație, rândul z al unei matrice pătrate este întotdeauna maparea z-ro a vectorului unitar (aici z = , ) Este foarte posibil să scrieți sub formă de matrice sistemul de ecuații ( ) ( , ) [X' y'] = [x y + [x-x y-y ] ( ) Totuși, prima parte a acestei ecuații nu este un produs pur matriceal În situațiile mai complexe în care rotația este combinată cu alte transformări, ar fi mai convenabil au un singur produs matriceal pentru fiecare transformare elementară La prima vedere, acest lucru pare imposibil dacă transformarea implică o operațiune de transfer Dar, așa cum vom vedea mai jos, cu ajutorul unei matrice de transformare * , acest lucru este foarte posibil Să începem cu un simplu transfer Fie punctul P(x, y) să fie transferat în punctul P'(x', y'), unde X' = X + O" [x' y' ] = [x y ] -sin sin^> -sinp cos C c k Cj = Xo - X COS Încheiem această secțiune cu trei observații: Fereastra nu trebuie neapărat să acopere întregul obiect ca întreg Dacă nu acoperă întregul obiect, atunci părțile obiectului care se află în afara ferestrei nu sunt desenate - trebuie tăiate Această operație se numește tăiere și este discutată mai detaliat în Secțiunea În cazul general, coeficienții /x și / sunt diferiți Pentru o diagramă cu bare, acesta este exact ceea ce aveți nevoie Dar nu este deloc potrivit atunci când rapoartele unghiulare de pe imagine trebuie să fie exact aceleași ca pe obiect În acest caz, cel mai mic dintre f și / ar trebui să fie ales ca factor de scalare Prin urmare, este recomandat O U LINII DE TĂIERE înlocuiți expresia ( ) cu formule bazate pe recalcularea coordonatelor relativ la centrele ferestrei și zona de ieșire Acest lucru va fi descris în paragraful Dimensiunile și poziția ferestrei nu sunt întotdeauna cunoscute dinainte Secțiunea va arăta cum sunt calculate în loc să fie setate de utilizator LINII DE TĂIERE Să presupunem că coordonatele lumii și ale ecranului sunt aceleași, adică fereastra și zona de ieșire sunt aceleași Prin urmare, în această secțiune, termenul "fereastră" poate fi înlocuit peste tot cu termenul "zonă de ieșire" Cu toate acestea, de obicei se consideră că tăierea se realizează pe marginile ferestrei, și nu pe zona de ieșire, iar acest lucru este luat în considerare în acest caz Pe fig arată dreptunghiul ABCD, care este o fereastră Toate segmentele de linie dreaptă vizibile trebuie să se afle în interiorul ferestrei Adică, atunci când desenați un segment de linie dreaptă, părțile acestuia aflate în afara ferestrei trebuie tăiate Procesul de tăiere ar trebui să fie automat Comenzi pentru desenarea unui triunghi PQR în fig sunt interpretate ca comenzi pentru a desena segmente de linii drepte Р'Р, PQ, QQ' Primul este desenat dreptunghiul ABCD, deci în loc de triunghi PQR, ar trebui trasată o polilinie P'PQQ' Orez Triunghi în curs de operație de tăiere Capitolul ALGORITMI D Deoarece sunt date doar trei puncte P, Q și R, perechea de coordonate de numere (xp>,yp,) trebuie calculată din valorile coordonatelor punctelor (xp, yp) și (xR, yR) Din fig Figura arată că panta segmentului PR poate fi calculată în două moduri, ceea ce duce la următoarea ecuație: Yr,-yr yR-yP Xp, - Xp XR - Xp Combinând această ecuație cu relația Yr'^max primim ^E-^pX 'max- 'r) Xp, = Xp + - p P XR~>'p De aici, coordonatele punctului P' sunt ușor de calculat dacă se știe că punctul final P se află în interiorul ferestrei, iar celălalt punct final R(Xp, yR) satisface inegalitățile xmin ^max Cu toate acestea, sunt mult mai multe situații de luat în considerare Varietatea mare de operații logice care trebuie efectuate pentru a rezolva această problemă fac ca problema tăierii liniilor să fie foarte interesantă din punct de vedere algoritmic Din fig este evident că nu este suficient să tăiați segmentul dreptei PQ în raport cu linia dreaptă CD Cohen și Sutherland au dezvoltat un algoritm pentru tăierea segmentelor de linie dreaptă, care este descris în Pascal în Newman și Sproul ( ) Aici, acest algoritm va fi prezentat în limbajul C Cu orice punct P(x, y) vom asocia un cod pe patru biți b b\ bc unde bі poate fi fie , fie (i = , , , ) Acest cod conține informații utile despre poziția punctului P în raport cu fereastra ABCD În limbajul C, expresiile condiționale produc valori Deci valoarea expresiei x xmax) /* punctul P la dreapta lui BC */ N = (y ymax) /* punctul P HadCD */ De fapt, doar nouă dintre cele combinații posibile de biți pot exista și sunt prezentate în Fig În limbajul C, valorile codului sunt generate de funcție int cod^(x, y) float x, y; { return (x xmax)" I (y ymax); } Pentru a înțelege această expresie, trebuie să știți că expresia b " n înseamnă că valoarea bitului b este deplasată cu n D cu oіЪo A ѳ oooh Orez Semnificațiile codului Capitolul ALGORITMI BIDIMENSIONALI poziții la stânga Pe lângă operatorul de deplasare pe biți ", există și un operator SAU pe biți, scris ca o bară verticală (I) Acest operator nu trebuie confundat cu operatorul logic SAU, notat cu bara verticală dublă (II) O modalitate mai puțin eficientă de a scrie expresia de mai sus într-o instrucțiune return ar fi (x x)* +(y y) v min max' -chtp' L Ltah' Această expresie este dată aici numai în scopul de a facilita percepția expresiei anterioare Funcția de cod descrisă aici va fi utilizată în funcția de clipare, a cărei sarcină este să analizeze un anumit segment de linie dreaptă și să deseneze numai acea parte a acestuia care este închisă în fereastra ABCD, dacă o astfel de parte există Această funcție funcționează după cum urmează Dacă cel puțin unul dintre codurile pentru punctele Pj și P conține un bit de , atunci fie Pp, fie P se mută din regiunea din afara ferestrei la una dintre limitele ferestrei sau la continuarea acesteia În acest din urmă caz, punctul va fi în continuare în afara ferestrei și va fi necesară încă o mișcare De exemplu, în cazul prezentat în fig , deplasarea din punctul P în punctul R trebuie urmată de trecerea din punctul R în punctul S Atunci poate fi necesară tăierea celuilalt capăt al segmentului, după cum se poate observa în fig Astfel, procesul de tăiere poate fi în mai multe etape, la fiecare pas distanța dintre punctele Pj și P scade Procesul se încheie imediat ce ambele puncte se află în fereastră Restul segmentului P^ va fi desenat Totuși, există un alt caz important în care bucla trebuie finalizată, și anume, când ambele puncte Pj și P sunt în afara ferestrei și în același timp pe aceeași parte a ferestrei Această situație nu poate fi distinsă la început, dar poate apărea în timpul procesului de tăiere Dacă punctele terminale ale unui segment de linie dreaptă sunt în afara ferestrei, atunci segmentul poate intersecta fereastra, dar poate fi și complet în afara ferestrei, așa cum se arată în Fig și Pe fig , puncte Pj și P nu sunt în același timp sub fereastră Dar în procesul de tăiere, punctele Q și S determină noile poziții ale punctelor Pj și, respectiv, P Deoarece acum ambele puncte sunt sub fereastră, putem concluziona că nu trebuie desenat nimic O astfel de decizie se ia pe baza analizei valorilor codeixl, yl) și codului (x , y ) Punctele Pt și P pe- LINII DE TĂIERE Orez Un segment de linie din afara ferestrei se deplasează pe aceeași parte a ferestrei dacă și numai dacă codurile lor conțin în aceeași poziție Pentru cele trei puncte A, Q și S, al treilea bit din stânga (M) din codurile lor este , în timp ce pentru punctul P , acest bit este Prin urmare, punctul P trebuie mutat în punctul S Mutarea punctului Pj la punctul Q nu este necesar, dar această mișcare se face pentru că nu dăunează, iar implementarea algoritmului este mai simplă La fel ca operatorul SAU pe biți, notat cu bara verticală (I), menționat mai devreme, C are un operator AND pe biți, notat cu & Reamintim că versiunile logice ale unor astfel de operatori sunt notate cu II și && Operatorii pe biți (notați cu I și &) au ca rezultat secvențe de biți, în timp ce operațiile logice au ca rezultat valorile simple sau , indicând "adevărat" și, respectiv, "fals" Nu numai , ci orice altă valoare, alta decât , reprezintă "adevărat" dacă este folosit într-o operație logică Astfel, construcția ciclică în timp ce (сі Іс ) va efectua acțiunile indicate de punctele de suspensie ( ) atâta timp cât evaluarea operației pe biți el I c va produce o secvență de biți care conține în orice Capitolul ALGORITMI BIDIMENSIONALI bit, adică atât timp cât există cel puțin un bit care conține în argumentele c sau c Pe de altă parte, operatorul if (el & c ) return ; determină o ieșire directă din funcție dacă secvența de biți rezultată din operația c &c conține cel puțin un -bit, adică atunci când ambele valori ale lui ci și c conțin un -bit în aceeași poziție Textul complet al funcției de clipare se află la sfârșitul următorului program, care decupează pentagoane imbricate /* Programul CLIPDEMO */ /★ Demonstrația algoritmului Cohen & Sutherland IIne-clipplng */ tfinclude "math h" float xmin- , xmax- , ymln- , ymax- ; maln() { Int I; float r, pl, alfa, phIO, phl, xO, yO, x , y , x , y ; pi- , *atan( , ); alfa- , *pl/ , ; phiO-OO; x - , ; y - , ; initgr(); /* Desenați chenarele ferestrei */ /* Fereastra este acum desenată */ mutare(xmin, ymin); draw(xmax, ymin); draw(xmax, colțuri); draw(xmln, ymax); draw(xmin, ymin); /* Desenați corect în cadrul ferestrei/ /* pentagoane concentrice */ /* În măsura în care este permis de limitele ferestrei, */ /* Sunt desenate de pentagoane regulate concentrice */ pentru(r- , ; r xmax)" I (y ymax); } clip(x , y , x , y ) float x , y , x , y ; { int c -code(x , y ) codul c (x , y ); float dx, dy; în timp ce (c c ) { if(c &c ) return; dx-x -x ; dy-y -y ; Dacă (c ) { Dacă (xKxmln) { y +- dy*(xmln-x )/dx; x -xmln; } else If (x >xmax) { y - dy*(xmax-x )/dx; x -xmax; } else If (yKymln) { x - dx*(ymin-y )/dy; y -ymin; } else If (y >ymax) { x' - dx*(ymax-y )/dy; y -ymax; } codul c (x , y ); } altfel { if (x xmax) { y - dy*(xmax-x )/dx; x -xmax; } else If (y ymax) { x - dx*(ymax-y )/dy; y -ymax; } c -code(x , y ); } } mutare(x , y ); draw(x , y ); } Pe fig arată rezultatul acestui program Orez Rezultatul programului CLIPDEMO Capitolul ALGORITMI BIDIMENSIONALI DIMENSIONARE AUTOMATĂ ŞI POZIȚII Pentru ca imaginea să fie desenată în limitele zonei de ieșire, trebuie mai întâi să efectuați tăierea pe fereastra specificată, așa cum este descris în paragraful , și apoi să reflectați fereastra și conținutul acesteia în zona de ieșire, așa cum se arată în Fig Pentru majoritatea aplicațiilor, această procedură este suficientă Dar acest paragraf va descrie o abordare ușor diferită, care diferă în următoarele aspecte: ) obiectul va fi desenat ca un întreg, deci nu este necesară tăierea; ) fereastra va fi determinată prin calcul și nu predeterminată; ) la oglindirea ferestrei pe zona de ieșire, același factor de scalare va fi utilizat pe ambele axe în direcția orizontală și verticală Din punctul rezultă că obiectul trebuie să fie finit Pentru majoritatea aplicațiilor, această limitare nu este severă, dar exclude panning Elementul poate fi implementat prin vizualizarea dublă a datelor care descriu obiectul În timpul primei scanări, limitele ferestrei sunt determinate xmin' Xmax' ^min" ^max* Desenarea se realizează în timpul celei de-a doua scanări În acest scop, vom folosi un fișier de pe disc pentru a elimina eventualele dificultăți din cauza limitării memoriei disponibile Punctul presupune că orice triunghi din imagine va fi similar cu triunghiul original de pe obiect, ceea ce înseamnă că relația unghiulară nu se va schimba atunci când este afișată Să presupunem că în fig avem un triunghi dat PQR, ale cărui coordonate ale vârfurilor sunt definite în sistemul de coordonate mondial Xp= , Xq= , xr= , J'p = , >'q = , yR=ll Atunci valorile calculate ale parametrilor ferestrei vor fi " - X = , minim maxim U = , la " , •pip Jmax Rețineți că punctele extreme ale obiectului sunt situate pe marginile ferestrei, ceea ce nu a fost cazul avut în vedere la paragraful DIMENSIONARE ȘI POZIȚIE AUTOMATĂ Pentru a oferi spațiu liber pe toate părțile ecranului sau ale foii de hârtie, trebuie să setați dimensiunile zonei de ieșire să fie puțin mai mici decât ar putea fi în realitate De exemplu, puteți seta xmin " , în loc de , Zona de ieșire va fi complet definită atunci când este selectată min = , Utyp = , x max Ymax = , = , Ca și în paragraful , calculăm factorii de scalare: ѵ - x max min xmav - x max mm , - , , - , = U-Utip , - , e max mm y "v -V , - , Ltah Lpt Ca factor de scalare f, alegem cea mai mică dintre valorile / x și / Amintiți-vă că toate distanțele sunt înmulțite cu un factor de scalare, deci dacă este ales un factor de scalare mai mare decât /x sau /y, atunci o parte a imaginii va depăși cu siguranță zona de ieșire În acest exemplu, obținem /-/X ib Capitolul ALGORITMI D În mod clar, același factor de scalare va avea ca rezultat un triunghi care este exact egal în lățime cu dimensiunea zonei de ieșire pe axa x, dar va exista spațiu liber în direcția y Este de dorit să-l distribuiți în mod egal între părțile inferioare și superioare ale zonei de ieșire Acest lucru se poate face alegând poziția centrului Yc în locul valorii minime Ymin, ca în paragraful , pentru a calcula constanta c Valoarea este calculată într-un mod similar xs x ' (xmjn+ Xmax^ = * ( - + , ) = , >t = - % Lax > =o- = = US " = ° = Сі=хс-/^с = - х \u d " - s \u d Noi- / noi \u d ; - x , - , Acum, pentru orice punct al obiectului (x, y), poziția punctului afișat (X, Y) este calculată prin formule X \u d / x + c \u d x- , Y \u d / -y + c \u d y - , Să scriem un program care va desena o imagine pentru care o fereastră nu poate fi setată în prealabil, iar acest exemplu va evidenția câteva aspecte noi ale limbajului C Folosim numere aleatorii pentru a genera o curbă de formă și dimensiune imprevizibile, unde, ca de obicei, curba este aproximată printr-un număr mare de segmente de linie dreaptă Scalarea și poziționarea automată vă permit să rezolvați o problemă care este practic insolubilă în orice alt mod Pentru fiecare segment, valorile coordonatelor x și y vor fi generate și scrise într-un fișier de pe disc Mai exact, vom scrie așa-numitele structuri care conțin triple codul x y unde perechea (x, y) definește coordonatele punctului în care stiloul ar trebui să se miște, iar parametrul codului poate lua valoarea sau , ceea ce înseamnă că starea stiloului este "ridicat" sau, respectiv, "coborât" Cu alte cuvinte, x y înseamnă toxe(X, Y) x y înseamnă draw(X, Y) unde X și Y sunt coordonatele ecranului corespunzătoare lumii DIMENSIONARE ȘI POZIȚIONARE AUTO coordonatele x și y În acest exemplu, vom folosi un program special care generează o curbă și scrie triplete de numere în fișierul A SCRATCH Execuția acestui program ar trebui să fie urmată de execuția programului general de trasare GEN PLOT, care citește triplete de numere de două ori Prima dată este să determinați parametrii ferestrei xtype, xmax, ymax, iar a doua oară este să efectuați efectiv operațiile de mutare și desenare a stiloului folosind coordonatele ecranului X, Y obținute prin conversia coordonatelor x, y în sistemul de coordonate mondial În program, vom începe să generăm curba din punctul de origine și vom deplasa de fiecare dată cu o unitate de distanță Există întotdeauna o direcție de curent^ și un unghi de rotație curent a În starea inițială, ambele aceste valori sunt egale cu zero Înainte de fiecare pas, unghiul a este mărit cu o valoare a unghiului selectată aleatoriu între - ° și + ° (întreg) Noul unghi de rotație rezultat este adăugat la direcția curentă^ pentru a obține noua direcție Să limităm valoarea maximă a curburii și să reducem șansa de degenerare într-un cerc Pentru a face acest lucru, modificăm algoritmul descris și setăm valoarea pentru unghiul a de fiecare dată când valoarea lui absolută depășește ° Apoi obținem următorul program: /* CURVGEN: generează o curbă aleatoare */ /* Generarea unei curbe aleatorii */ #include "math h" #include "stdlo h" #include "tlme h" /* Acest fișier antet este necesar când se apelează */ /* funcția de timp pentru a obține un număr aleator */ /* Acest fișier antet definește macro-ul abs */ maln() {Int i N- ; float x ) O, yO O, xO, yO, phl, direction(); pfopen(); mutare(x, y); pentru ( - ; KN; i++) {xOx; yO-y; direcția phi(); x-xO+cos(phi); y-yO+sin(phi); draw(x, y); } pfclose(); } - Capitolul ALGORITMI BIDIMENSIONALI float dlrection() { static int phl-O, alfa-O, primul- ; /♦Variabilele statice sunt inițializate numai în primul caii! */ /*Variabilele statice sunt inițializate doar la primul apel!*/ float pl- ; semințe lungi Int; Dacă (primul) { primulH); tlme(&seed); srand((lnt)seed); } alfa-*-rand( )% - ; dacă (abs(alfa)> ) alfa-O; phl-H-alfa; return ((float)phl*pi/ ); } FIȘIER *fp; struct {float xx; plutește yy; intII;}s; pfopen(){ fp-fopen("a scratch", "wb"); } /* wb; netradus */ /* wb - fișier deschis pentru scriere în modul binar */ ptoѵe(x, y) float x, y; {s xx-x; s aa-a; s ll-O; /* - rep up */ /* - pen up */ fwrlte(&s, dimensiunea s, , fp); } pdraw(x, y) float x, y; {s xx-x; s aa-a; s il- ; /* - stilou jos */ /* - stilou jos */ fwrlte(&s, dimensiunea s, , fp); } pfclose() { fclose(fp); } Funcția de direcție arată un apel la funcția de inițializare a generatorului de numere aleatoare srand Argumentul său de semințe specifică valoarea de semințe pentru generarea de numere aleatorii Funcția de timp este utilizată pentru a atribui o valoare argumentului semințe pe baza orei curente În acest fel, vor fi generate diferite curbe de fiecare dată când programul este repornit Funcția rand produce un întreg mare nenegativ Este convertit într-un număr întreg în intervalul de la la utilizând restul împărțirii cu , adică în program, care "include" fișierul antetelor funcției standard de intrare/ieșire Limbajul C face distincția între I/O formatat și neformatat Cele mai frecvent utilizate funcții sunt: DIMENSIONARE ȘI POZIȚIONARE AUTO scanf: intrare terminal formatat printf: ieșire formatată către terminal fscanf: intrare formatată dintr-un fișier fprintf: ieșire formatată în fișier fread: intrare neformatată dintr-un fișier fwrite: ieșire neformatată într-un fișier I/O formatat se ocupă de caractere care pot fi citite: există o structură de linii similară liniilor de pe o pagină tipărită Datele neformatate au aceeași structură în care sunt stocate în memorie De exemplu, numerele întregi sunt reprezentate printr-un număr fix de biți În cazul nostru, I/O brut a fost folosit pentru a îmbunătăți eficiența Fișierul este "deschis" prin apelarea funcției fopen Al doilea argument poate fi fie șirul "r" pentru a iniția o citire, fie șirul "u>" pentru a scrie (Pentru I/O brute, unele compilatoare necesită șiruri altele decât "r" sau "w" ) Variabila $ denotă o structură care conține trei numere, x, y și cod, care vor fi scrise în fișier "Indicatorul" &s la aceasta variabila este primul argument al functiei fwrite Intrarea poate fi considerată ca indicând adresa variabilei $ Al doilea argument, sizeof s, este egal cu dimensiunea unei structuri $, al treilea argument, egal cu , înseamnă numărul de structuri care trebuie scrise Al patrulea argument, /p, este indicatorul fișierului obținut din apelul la fopen Acum putem lua în considerare un program general care va citi tripleți, va determina dimensiunile ferestrei, va converti coordonatele lumii în coordonatele ecranului și, în final, va desena o imagine într-o zonă de ieșire dată De asemenea, desenează colțuri mici în colțurile zonei de ieșire și un punct în mijlocul marginii de jos, astfel încât imaginea abstractă să poată fi orientată în raport cu partea de sus și de jos (dacă punctul și colțurile sunt nedorite, acestea pot fi ușor exclus) Rețineți că acest program deschide fișierul A SCRATCH de două ori Operațiunile de închidere și deschidere "derulează" fișierul la început /* GE MPLOT: */ /* Un program general de reglare și grafică */ /* Fișierul A SCRATCH conține datele de intrare */ ** Capitolul ALGORITMI BIDIMENSIONALI /* Program generic de dimensionare și desen */ /* Fișierul A SCRATCH conține date de intrare */ #include "stdio h" principal() { float X, y, xmin, xmax, ymin, colțuri, X, Y, Xmlri, Xmax, Ymin, Ymax, fx, fy, f, xC, yC, XC, YC, c , c ; FIȘIER *fp; struct {float xx; plutește yy; IntII;}-uri; fp-fopen("a scratch", "rb"); /* Parametrul "rb" este necesar pentru citirea în modul binar */ /* necesar în Lattlce C pentru modul binar */ xmln-ymin- e ; xmax-ymax-xmin; whlle (fread(&s, slzeof s, , fp)) {xs xx; da aa; Dacă (x xmax) xmax-x; dacă (y ymax) ymax-y; } fclose(fp); lnlt viewport(&Xmln, &Xmax, &Ymin, &Ymax); fx-(Xmax-Xmin)/(xmax-xmIn); fy-(Ymax-YmIn)/(ymax-ymin); f-(fx Uo> PlCxpjj) Pn - Y^ apoi partea curbei B-spline dintre două puncte consecutive P* iar P/+ se obține prin calcularea funcțiilor x( și y(t) pentru schimbarea t de la la x(t) = {(a t + ax} +a y(O = {(^ + * )/ + M+yo Aceste ecuații conțin următorii coeficienți: az = (-x^ + x{ - x/+ + / a = ^Xi- " d + Xi+ > / ( ) ai = C \u d P, + Xv * \u d x ( ) \u d (xA + xv + xc) / unde xb* denotă valoarea calculată a coordonatei x pentru punctul B Considerând același punct ca aparținând segmentului BC, obținem A = P/+p B-Pp C = P|+| xv * - x ( ) \u d (xd + xv + xs) / Capitolul ALGORITMI D Aceasta arată că ambele metode de calcul a valorii lui x dau același rezultat, ceea ce înseamnă continuitatea funcției x (t) în punctul B Diferențiând x ( de două ori, găsim derivatele x' ( și x "( ) Înlocuind în ele valorile t- și te , așa cum sa făcut pentru x(t), se va putea verifica dacă derivatele sunt continue în punctul B Deoarece funcția y(t) și primele ei două derivatele sunt, de asemenea, continue, devine clar că curba B-spline este foarte netedă Pentru a calcula orice segment de curbă între punctele Pf și Pzh, se folosesc și punctele Pm și P^ Rezultă de aici că primul segment al curbei va fi situat între punctele Pj și P , iar ultimul segment între punctele Pj și P p, astfel încât punctele de început și de sfârșit ale întregii curbe să fie situate în apropierea Pj și Pp p dar nu lângă Po și RL Programul de mai jos citește numere P * * Wo Y} x "K n*n din fișierul CURV DAT În ieșire, fiecare dintre cele n + puncte este notat cu un marker sub forma unei cruci Apoi este trasată o curbă B-spline /*CURVFIT: netezirea curbei folosind B-spline */ /* Ajustarea curbei folosind B spllnes */ #include "stdio h" #definiți MAX tfdeflne N principal() { float x[MAX], y[MAX], eps- , , XY t xA xB xC xD yA yB, yC, yD, aO, a , a , a , b , s, b , b ; int n, i, j, flrst; FIȘIER *fp; fp-fopen("curv dat", "r"); if (fp NULL){ printf /* "Nu există niciun fișier curv dat" */ ("Fără fișier curv datXn"); exlt( ); } fscanf(fp, "%d", &n); formă); i y = y + ((Ni)R/N) sin^> Luați valori fixe y , R, N', de exemplu, = , y = , , R = , N= Fie că i variază în intervalul de la la N - Fie pentru fiecare valoare a lui i, unghiul ia succesiv valorile °, °, °, , ° Scrieți un program pentru a desena un set de de triunghiuri Vârfurile primului triunghi sunt date în punctele ( , ), ( , ), ( , ) Fiecare triunghi următor se obține prin rotirea celui precedent cu un unghi de ° în jurul punctului ( , ) Implementați un algoritm recursiv de tăiere a liniilor bazat pe bisectarea unui segment de linie Să presupunem că sunt date o fereastră și un segment de linie dreaptă PQ Dacă ambele puncte P și Q sunt în interiorul ferestrei, atunci segmentul este desenat complet Această condiție acoperă și cazul în care un punct se află în interiorul ferestrei, iar al doilea este la graniță Există situații în care puteți determina cu ușurință că nu trebuie desenat nimic Execută-le pe cont propriu În toate celelalte cazuri, este necesar să găsiți punctul M în mijlocul segmentului PQ și să aplicați recursiv această procedură la segmentele PM și MQ Asigurați rezistență la marginile ferestrelor și investigați efectul numărului de bisecții ale segmentelor de linie asupra acestei stabilități Proiectați și implementați un algoritm pentru tăierea părților de linii dacă fereastra este dată în formă de triunghi Modificați programul STELE din Secțiunea , astfel încât stelele să nu se suprapună Utilizați numere aleatoare pentru a genera o secvență de, să zicem, de puncte și aplicați un program de netezire pentru a desena o curbă prin aproximativ acele puncte Scrieți un program care generează un arbore mai realist decât cel prezentat în Fig Un exemplu de imagine este prezentat în fig Orez Copac capitolul INSTRUMENT GEOMETRIC PENTRU ALGORITMI GRAFICI D VECTORI Este necesar să aveți o anumită pregătire matematică pentru a înțelege și a scrie programe pentru grafică tridimensională Această carte nu se dorește a fi un manual de matematică și se presupune că cititorul este deja familiarizat cu tehnicile matematice folosite în acest capitol, în special cu vectorii și determinanții Un vector este un segment direcționat al unei linii drepte, caracterizat doar prin lungimea și direcția sa Pe fig prezintă două reprezentări ale aceluiași vector a = PQ = b = RS Prin urmare, vectorul nu se modifică în timpul translației Pe fig punctul de început al vectorului b este același cu punctul final al vectorului A Apoi suma vectorilor a și b este definită ca vectorul c tras de la punctul de pornire al vectorului a până la punctul final al vectorului b, ca să putem scrie c = a + b Orez Vectori egali VECTORI Lungimea vectorului a se notează cu I a I și este egală cu distanța dintre punctele sale de început și de sfârșit Un vector cu lungime zero se numește vector zero și se notează Notația -a este folosită pentru un vector de lungime I a I, a cărui direcție este opusă direcției vectorului a Pentru orice vector a și un număr real c, vectorul ca are lungimea I c a I Dacă a = sau c = , atunci ca = , în caz contrar vectorul ca coincide în direcția cu vectorul a dacă c > , sau are direcția opusă dacă c ij = ji = jk = kj = ki = ik = O Echivalând a = b în ecuația ( ), obținem a-a = I a , astfel încât I a I \u d V la-a I Proprietăți importante ale produsului punct: s (Li-ѵk \u d sL (și-v) (cu + Lv) • w = cu• w + Lv- w uv = vu u • u = O numai dacă u = O DETERMINANȚI Produsul scalar al vectorilor u = [w u u ] și v= [Vj v v ] poate fi calculat ca U'v = w v + w v + u v Acest lucru este ușor de demonstrat prin rescrierea părții drepte a ecuației u-v = (Wjі + zz j + zz k) • (gr + v j + v k) ca sumă a nouă produse scalare și analizându-le pe baza expresiei ( ) DETERMINANTE Înainte de a descrie produsul vectorial, să acordăm atenție determinanților Pentru a rezolva următorul sistem de două ecuații liniare: [a x + b y = c este necesar să înmulțim prima ecuație cu coeficientul b , iar a doua - cu coeficientul -bx și să adunăm, apoi obținem (ai ^ # ^ % ^ ^ ^ După aceea, puteți înmulți prima ecuație cu -a , iar a doua - cu topor și o puteți adăuga și ea Drept urmare, obținem {axb -arbx)y=axc -a cx Dacă axb - a bx nu este egal cu zero, atunci puteți efectua împărțirea și găsirea b'fC-t blC*) a^c^ a^c i > y= -G~G ( , ) (axb ~a bx) J (axb -a bx) Expresia din divizor se poate scrie sub forma topor bx a b În acest caz, se numește determinant de ordinul doi Prin urmare, a * a b a\^ a ^\ CU folosind determinanți, se poate scrie ecuația ( ) sano în formă D\D Capitolul INSTRUMENT GEOMETRIC Unde £> = J * a b ^ = C L C b D = fll C a c Rețineți că D se obține prin înlocuirea coloanei /-a din D cu partea dreaptă a sistemului de ecuații ( ) (/ = sau ) Acest mod de rezolvare a unui sistem de ecuații liniare se numește "regula lui Cramer" Această metodă este potrivită nu numai pentru un sistem de două ecuații (deși din punct de vedere al timpului mașinii se dovedește a fi foarte costisitoare pentru sistemele mari) Definim determinantul de ordinul trei sub forma unei ecuații D= a\ b\ C a ^ C a ^ C = J b c bz c ~а b\ bz c\c + a bx CX b c iar determinantul de ordinul al patrulea D= d * C d\ a ^ C C a ^ c d a b C d c d bx cx dx C dl Cj d{ = a ^ c d ~a ^ c d +az~ ^ c d ~a ^ c d b c d b c d b c d ^ c d și așa mai departe Determinanții au multe proprietăți interesante, unele dintre care sunt enumerate mai jos Valoarea determinantului nu se va modifica dacă rândurile sunt scrise ca coloane în aceeași ordine, de exemplu: aX a b ah ii a b Dacă înlocuim reciproc două rânduri (sau două coloane), atunci valoarea determinantului va fi înmulțită cu - : J * C a ^ C a ^ C J ^ C a ^ C a ^ C Dacă orice rând (sau coloană) este înmulțit cu un factor, atunci și valoarea determinantului va fi înmulțită cu acel factor De exemplu: DETERMINANTE sah cbx a b ' va fi descris de către ecuaţie X y Z XI Y! Z, X y Z xs Uz z " Capitolul INSTRUMENT GEOMETRIC VECTOR ART Produsul încrucișat al doi vectori a și b este notat a*b și este egal cu vectorul ѵ, care are următoarele proprietăți Dacă a = cb pentru un scalar c, atunci v = În caz contrar, lungimea vectorului v este I ѵ I \u d I a I ІЫ sin y unde y este unghiul dintre vectorii a și b, iar direcția vectorului v este perpendiculară pe ambii vectori a și b și este astfel încât a, b, v, în această ordine, formează un triplu dreptaci Aceasta din urmă înseamnă că dacă vectorul a se rotește printr-un unghi y j)k care poate fi scris ca a a i + a a i + " a & ^ b b{ ^ sau într-o formă mai ușor de reținut: a x b = ij k a\ a a ^ ^ VECTOR ART Orez Produsul încrucișat al lui k-a x b Aceasta este mai mult o notație mnemonică decât un determinant real, deoarece elementele primului rând sunt vectori, nu numere Dacă vectorii a și b indică laturile adiacente ale paralelogramului, ca în fig , atunci aria acestui paralelogram este egală cu lungimea vectorului a x b Aceasta rezultă direct din notația I a x b I \u d I a b I sin y Pe fig vectorii a și b se află în planul care trece prin axele x și y Să presupunem că axa z iese din foaia de hârtie la Capitolul INSTRUMENT GEOMETRIC cititor și corespunde sistemului de coordonate corect Apoi pentru spațiul tridimensional ae [aj a ], fi [Aj b ] ijk a\ a b o a*b = a ' b Astfel, vectorul a x b va avea aceeași direcție ca și vectorul k, dar numai dacă determinantul are un semn pozitiv Aceasta impune condiția ca vectorul a, rotit spre vectorul b cu un unghi mai mic de °, să se rotească în direcția pozitivă (în sens invers acelor de ceasornic) dacă și numai dacă D > Vom folosi această metodă mai jos pentru a determina dacă vârfurile triunghiului A , B, C în sens invers acelor de ceasornic atunci când sunt listate în această ordine Pe fig avem u=[zz u ] = AB, ve [Vj ѵ ]vAC Orez Punctele A, B, C sunt ocolite în sens invers acelor de ceasornic DESCOMPUNEREA POLIGONILOR ÎN TRIANGURI HA W A xv xs Noi D= HA Marea Britanie hv hl Uv~Uk xs "xA noi~ya ХВ~ХА >'v~>'a xs~ha Us~Uk I V ȘI V Vârfurile A, B, C, în această ordine, sunt ocolite în sens invers acelor de ceasornic dacă, și numai dacă, vectorul și se rotește în direcția vectorului v cu un unghi y - punctele A, B, C sunt ocolite în sens invers acelor de ceasornic; D •••' Po^n- POATE se extinde dincolo de poligon Acum să creăm un program care să citească coordonatele vârfurilor poligonului și să împărțim poligonul în triunghiuri Este necesar ca vârfurile să fie specificate în sens invers acelor de ceasornic De exemplu, la poligonul din fig ( ) secvența P P P P P P P va fi corectă, iar secvența P^P^ P^Po va fi nepotrivită Pentru un poligon cu n vârfuri, mai întâi este specificat numărul de vârfuri, iar apoi n perechi de coordonate ale tuturor nodurilor sunt listate secvenţial în ordinea inversă a acelor de ceasornic în jurul poligonului Ca urmare, se va obține un desen de poligon cu diagonale desenate care împart complet întregul poligon în triunghiuri Înainte de a desena diagonalele, trebuie să vă asigurați că toate diagonalele se află complet în interiorul poligonului Să presupunem că Pb , P-, P/+ denotă trei vârfuri învecinate și vom presupune că P t = Pn- și Pn = Po, astfel încât să putem lua în considerare cazurile în care / = și i = n - Reamintim că vârfurile trebuie listate în ordinea de parcurgere în sens invers acelor de ceasornic În acest caz, P va fi un vârf convex dacă și numai dacă trei vârfuri P ^ p P-, P + în această ordine vor ocoli în sens invers acelor de ceasornic În acest scop, folosim o metodă care funcționează bine pentru o clasă largă de poligoane Cu toate acestea, pot exista poligoane pentru care această metodă dă o eroare Vezi, de exemplu, exercițiul de la sfârșitul acestui capitol Capitolul INSTRUMENT GEOMETRIC săgeți care urlă Ca contraexemplu, luăm în considerare Fig , în care parcurgerea triplei * * ^ se realizează în sensul acelor de ceasornic Vârful P este neconvex, iar diagonala PjP se află în afara poligonului Astfel, diagonala P^ P^ poate fi candidată pentru scopurile noastre numai dacă trei vârfuri P/ y^), P / Xp y ), P / + (x + , y / + ), în această ordine, sunt ocolite în sens invers acelor de ceasornic, adică dacă D= хі уі Vi Vi > Această condiție este necesară, dar, din păcate, nu este suficientă, așa cum se arată în Fig Aici punctele Po, Pp P sunt ocolite în această ordine, în sens invers acelor de ceasornic Orez Diagonala Р Р este parțial în afara poligonului DESCOMPUNEREA POLIGONILOR ÎN TRIANGURI săgeți, dar segmentul P P nu poate fi folosit pentru a împărți un poligon în triunghiuri Această situație poate fi evitată dacă se ține cont și de lungimea diagonalelor Vom alege cea mai scurtă diagonală P^P^ pe care o poate avea vârful convex Pf între punctele P^ și P^ Această diagonală este folosită pentru a tăia triunghiul P Pf-P^ • Apoi poligonul rămas se verifică în același mod Po, Pp , Pf p P/+p , P^! și așa mai departe Din punct de vedere tehnic, acest lucru este implementat prin introducerea unui tablou întreg ѵ , , ѵ р care conține numărul de vârfuri ale poligonului rămas Mai întâi, setăm m=n și ѵ =i (z= , , , m- ) De fiecare dată când triunghiul este tăiat, numărul m scade cu unu Dacă un astfel de program este destinat utilizării practice, atunci este de dorit să se efectueze o serie de teste pentru validitatea datelor de intrare Se poate spune că un program care respinge în mod fiabil orice set de date inutilizabil este robust În situația noastră, este necesar să testăm dacă setul dat de coordonate de puncte este într-adevăr descrie în general poligonul De exemplu, secvența întrucât ocolirea punctelor în ordinea dată va duce la situația prezentată în Fig , dar o astfel de cifră nu poate fi acceptată ca poligon Orez Rezultatul secvenței nevalid Capitolul INSTRUMENT GEOMETRIC Alte verificări de luat în considerare includ: - numărul maxim de puncte și, de exemplu, n -NMAX) { printf /* "n prea mare" */ ("numărul n este prea mare"); ieșire( ); } pentru H); i ) { min dlag-BIG; DESCOMPUNEREA POLIGONILOR ÎN TRIANGURI pentru (Yu; Km; I++) { h- (I О ? m- ; i- ); j- (I t- ? O: + ); Dacă (counter clock(h, I, j, &diag) && diag e- ; } draw poligon() { Int I; mutare(x[n- ], y[n- ]); pentru (i- ; Kn; I++) draw(x[l], y[i]); } liniuță(x , y , x , y ) float x , y x , y ; { Int I, k; float xdif-x -x , ydif-y -y , pitchO-O , dx, dy; k - * (int)celulă(sqrt(xdlf*xdlf+ydlf*ydlf)/pltchO)+ ; dx-xdif/k; dy-ydlf/k; pentru (I- ; Kk; I+- ) { mutare(x +i*dx, y +l*dy); draw(x +(l+ )*dx, y +(i+ )*dy); } } Următoarea secvență de introducere - Capitolul INSTRUMENT GEOMETRIC Orez Rezultatul programului POLI TRIA va produce imaginea prezentată în Fig COORDONATE OMogene Această secțiune tratează problema interesantă din punct de vedere matematic a proiecției în perspectivă și, prin urmare, este adesea inclusă în cărțile despre grafica pe computer D Numai proprietățile geometrice vor fi discutate aici Deși argumentele sunt relativ simple, ele sunt destul de lungi și teoretice Scopul este de a înțelege mai bine notații precum [x y ] și [x y z ] Cititorii interesați de aspectele practice ale programării grafice pot sări peste această secțiune fără consecințe grave pentru înțelegerea restului cărții Puteți reveni la studiul acestui paragraf într-o etapă ulterioară Secțiunea a folosit deja notația [x y ] pentru a desemna o matrice constând dintr-un singur rând, uneori numit un vector rând O astfel de înregistrare poate fi considerată și ca un caz special al înregistrării [xy w], unde numerele x, y, w se numesc coordonate omogene Aceste trei numere de coordonate omogene sunt folosite pentru a desemna un punct din spațiul bidimensional În geometria proiectivă, coordonatele omogene au fost folosite cu mult înainte de apariția și răspândirea graficii pe computer În capitolul vom lua în considerare transformările de perspectivă De fapt, ele sunt proiecții centrale, explorate în diverse cărți despre geometria proiectivă Aici discutăm doar pe scurt câteva probleme din această problemă fascinantă, dar destul de greu de înțeles COORDONATE OMogene ramură a matematicii, evitând definițiile formale și dovezile riguroase ale teoriei Pe fig avem axa x și axa w, astfel încât punctul este dat de o pereche de coordonate (x, uO Orice punct P (x, y) care nu se află pe axa x are proiecția sa centrală P' (X, ), definit ca punctul de intersecție linia dreaptă OP cu o dreaptă Z descrisă de ecuația w = Punctul de origine este centrul proiecției Segmentul de dreaptă PO poate fi considerat ca un fascicul de lumină de la obiectul P la ochiul situat în punctul O Pentru punctele Q( , u și Q '( , ) obținem două triunghiuri similare OPQ și OP'Q', astfel încât X P'Q' PQ x OQ' OQ w Toate punctele (x, uO cu proprietatea x = wX se află pe linia OP și au aceeași proiecție P' Dacă ne interesează doar proiecțiile pe dreapta Z, și nu valorile reale ale lui x și iv, atunci contează doar raportul lor Este destul de natural să folosiți o singură coordonată X în loc de o pereche (X, ), dacă sunt luate în considerare doar punctele de pe o linie dreaptă Z Dacă tot trebuie să utilizați o pereche de coordonate , atunci orice pereche de numere (x, id, care satisface condiția X, dacă acceptăm un astfel de acord În sens geometric, perechea de coordonate (uX, u>) a oricărui punct P, altul decât O, de pe linie Orez proiecție centrală D la naiba Capitolul INSTRUMENT GEOMETRIC linia OP' poate servi ca desemnare a punctului P' Acest lucru se realizează în cazul utilizării coordonatelor omogene În general, orice punct (Xp X , în spațiul n-dimensional este scris ca punct (wXp wX , uXn, w) în (n+D-spațiul dimensional, unde w este orice număr real diferit de zero Acest grup de m Numerele + definesc coordonatele omogene ale punctului original din spațiul n-dimensional Cititorul poate fi familiarizat cu definiția unei proiecții, care este o mapare a unui set de puncte din spațiul (m + )-dimensional într-unul în m -spațiu dimensional de la spațiul n-dimensional la spațiul (m + )-dimensional Punctul ( , ) din spațiul bidimensional, de exemplu, poate fi scris în coordonate omogene ca ( , , ) sau ca ( , , ) și așa mai departe Deși toate operațiile se referă la spațiu bidimensional, aceste triple pot fi luate ca puncte în spațiul tridimensional Evident, notația (X, Y) este notația obișnuită pentru un punct din spațiul bidimensional pentru un punct (X, Y, ) din spațiul tridimensional Punctul P' este central pentru proiecție pentru orice punct P (x, y, u>), dacă x / w \u d X și y / w \u d Y Ca și mai înainte, punctul de origine O va fi centrul proiecțiilor și toate punctele sunt proiectate pe avionul w \u d Pentru a defini termenul omogen (coordonate), folosim ecuația aX+SY+c= ( , ) descriind o linie dreaptă în spațiul bidimensional Înlocuind X și Y cu x/w și y/iv, obținem a(x/w) + b(y/w) + c = sau ax + sy + civ = ( , ) Ecuația ( ) este de obicei numită omogenă, deoarece are aceeași structură în ceea ce privește ax, by, cw Prin urmare, numerele x, y, w se numesc în mod natural coordonatele omogene ale punctului (X, Y) Dacă presupunem din nou că spațiul bidimensional este situat în planul w= în sistemul de coordonate xyw, atunci ecuația ( ) descrie planul care trece prin origine și linia dreaptă dată COORDONATE OMogene Dacă presupunem că notația (x, y, u>) este folosită ca o altă formă de notație pentru (x/ u>, y/u>), atunci va fi necesar să cerem ca valoarea lui w să nu fie egală cu zero Cu toate acestea, în acest caz, coordonatele omogene nu au avut avantaje față de coordonatele obișnuite, iar cititorul își putea exprima în general îndoielile cu privire la prezența oricăror avantaje De exemplu, luați în considerare un sistem de două ecuații liniare, fiecare dintre acestea descriind o linie dreaptă în spațiu bidimensional: a X + b y + c ^ Dacă două drepte sunt paralele, atunci ele nu se intersectează și nu există o pereche de numere (X, Y) care să satisfacă sistemul ( ) Astfel, pentru a găsi un punct comun pentru linii drepte, este necesar să se aplice regula cu unele excepții, care nu este în întregime elegantă Când înlocuiți coordonatele X și Y cu coordonate omogene x, y, u, situația se va îmbunătăți oarecum |a x + A y + c u>= ( ) (a x + b y + c w = Ecuațiile din sistemul ( ) pot fi interpretate ca plane care trec prin punctul de origine O Acest sistem are cel puțin o soluție trivială x = y = w = Pentru interpretarea geometrică, setăm valori specifice pentru coeficienți Să înlocuim, de exemplu, sistemul ( ) cu sistemul X + ZU- = X + U- \u d O descriind două linii drepte paralele prezentate în fig Atunci sistemul de ecuații ( ) poate fi înlocuit cu G x + y- w= (Z Yu a) [ x + boo - u > = ( ) Acest sistem este echivalent cu sistemul x + Zy = w = astfel încât soluția să fie formată din toate tripletele ( L, - L, ), unde k este orice număr real În spațiul tridimensional, aceste puncte formează o linie dreaptă care trece prin punctele O și ( , - , ), Capitolul INSTRUMENT GEOMETRIC în plus, această dreaptă este linia de intersecție a celor două plane dată de ecuațiile ( ) Revenind la spațiul bidimensional al planului , reamintim că fiecare punct (X, Y) este asociat cu o dreaptă (u>X, wY, w) în spațiul tridimensional Pentru valorile diferite ale lui u, această asociere este aproape banală Acum va deveni clar un motiv foarte important pentru utilizarea coordonatelor omogene Pentru fiecare linie dreaptă din spațiul D, să adăugăm un obiect numit punctul de la infinit Acest punct la infinit nu poate fi notat în sistemul de coordonate obișnuit, dar în sistemul de coordonate omogene se poate De exemplu, un punct la infinit pe o dreaptă descrisă de ecuația ( IO a) este scris ca ( , - , ) sau ca orice triplă ( L, - A, ) pentru k diferit de zero Deoarece aceste triple sunt soluții ale sistemului de ecuații ( ), atunci punctul de la infinit poate fi considerat ca punctul de intersecție a două drepte paralele prezentate în Fig Un punct la infinit este considerat a fi în spațiu bidimensional După cum am văzut, fiecare punct din spațiul D este asociat cu o linie dreaptă din spațiul D, așa că este de dorit să aflăm cu ce linie dreaptă este asociat punctul de la infinit ( , - , ) Deoarece această linie trebuie să treacă prin punctul de origine ( , , ), atunci linia dorită trebuie să fie o linie care trece prin punctul O și punctul ( , - , ), adică situată în planul ѵv= Este destul de rezonabil să numim punctul ( , - , ) un punct infinit de distanță, deoarece poate fi considerat ca un punct limită ( , - , u>) deoarece u> tinde spre zero și acest triplu în coordonatele omogene este echivalentă cu punctul ( /w, - /u>, ), care este foarte îndepărtat pentru valori mici ale lui w Introducerea unui punct la infinit ne permite să afirmăm că oricare două drepte distincte se intersectează într-un punct În mod similar, în geometria proiectivă se poate argumenta că oricare două plane diferite au o linie de intersecție Dacă planurile sunt paralele, atunci toate punctele acestei linii de intersecție în coordonate omogene sunt scrise ca (x, y, z, ) Nu vom lua în considerare această problemă în detaliu, ci mai degrabă ne vom întoarce la spațiul bidimensional și vom arăta alte manifestări noi ale posibilităților oferite de coordonatele omogene COORDONATE OMogene L = În coordonatele obișnuite, neomogene, o transformare liniară în spațiul bidimensional poate fi scrisă ca [X' Y'] = [X Y] Deoarece [ ] А = [aj а ] şi [ ] Л = [й /> ], rândurile matricei A reprezintă, respectiv, punctele [ ] şi [ ] Indiferent de modul în care este definită matricea A, originea O nu se modifică, astfel încât [ ] A = [ ], deci operația de transfer nu poate fi exprimată în acest fel Totuși, în coordonate omogene, un punct din spațiul bidimensional este dat de un triplu (x, y, z) și transformarea se scrie ca y iv]L *h] bz S x [x' y w' ]i [x a C l = Avem a ^ c Orez linii drepte paralele Capitolul INSTRUMENT GEOMETRIC Deoarece punctul [ ] este un punct la infinit pe axa x, primul rând [Oj a a ] al matricei A este o mapare a acestui punct la infinit pe axa x În mod similar, al doilea rând al matricei b b ] va fi afișarea punctului la infinit pe axa y Deoarece [ ]А Hq с с ] atunci putem presupune că al treilea rând al matricei [q c c ] este o mapare a punctului de origine [ ] Aceasta înseamnă că coordonatele omogene permit exprimarea oricărei transformări prin înmulțirea matricei De fapt, acest lucru nu este nimic nou și această proprietate a fost deja folosită în Secțiunea Cu toate acestea, operația de mutare nu este singura proprietate nouă oferită de o astfel de multiplicare a matricei Poate fi folosit pentru a converti linii drepte paralele în linii care se intersectează Vom arăta acest lucru prin exemplul conversiei unui cadran dreptunghiular complet într-un triunghi Pe fig prezintă un triunghi arbitrar, ale cărui vârfuri sunt date într-un sistem de coordonate dreptunghiular de punctele A (ap a ), B (Z> p b ) și C (cp c ) Să adăugăm o a treia coordonată, egală cu , ca mijloc formal de formare a coordonatelor omogene Excepție fac punctele infinit îndepărtate ( , , ) și ( , , ) de pe axele de coordonate, pentru care a treia coordonată omogenă este egală cu zero Formăm o matrice I, astfel încât înmulțirea matricei [*' y' z' ] = [x y z]M va mapa punctul O la punctul C, punctul ( , , ) la punctul A și punctul ( , , ) la punctul B Pentru orice valori diferite de zero ale a,/î, punctele A, B și C vor avea puncte ale imaginii dorite a triunghiului, dacă L/ = aah aa a ?b p noi{ noi y ( , ) Acest lucru este ușor de verificat, deoarece avem relația evidentă [ ] M = [aa aa a] a cărui parte dreaptă este doar o altă desemnare pentru punctul A din fig La prima vedere, se pare că constantele a, /?, y din ecuația ( ) nu sunt necesare și pot fi stabilite COORDONATE OMogene 'Orez Cuadrant transformat în triunghi egale cu Cu toate acestea, acestea sunt necesare pentru a putea seta o mapare predeterminată U' pentru așa-numitul punct unitar U Punctul U' poate fi ales oriunde în triunghi Deoarece punctul U(l, , ) este mapat pe punctul U'(ttp u , m ), avem [ și ] este o notație prescurtată pentru următorul sistem de trei ecuații liniare în raport cu variabilele a, /?, y: a a + d ^ + c y = u a a + b/ + c y = u a + [i + y = Sistemul are o soluție unică, care rezultă din faptul că triplele (ap a , a ), (bv b , Z> ) și (cp c , c ) denotă vârfurile unui triunghi Rezolvarea acestui sistem pentru a, / , y și înlocuirea rezultatului în ecuația ( ) dă valoarea dorită a matricei A/ Liniile drepte sunt convertite în linii drepte, dar paralelismul lor nu este păstrat De exemplu, linia verticală care trece prin punctul U din fig se transformă într-o linie dreaptă care trece prin punctele B și U' Aceasta oferă un instrument geometric pentru găsirea proiecției P' pentru punctul P Având în vedere valoarea calculată a matricei A/, punctul de proiecție P' poate fi găsit analitic ca produs [ ]M Capitolul INSTRUMENT GEOMETRIC Rețineți că toate punctele de la infinit sunt mapate pe segmentul AB și toate punctele de la infinitate de linii paralele sunt mapate pe un singur punct de pe segmentul AB Acest lucru dă motive să considerăm întreg cadranul ca o perspectivă plată, în care triunghiul ABC este imaginea, iar segmentul AB este linia orizontului Aceasta completează analiza coordonatelor omogene Cititorul poate găsi alte proprietăți interesante ale coordonatelor omogene în manualele de geometrie proiectivă, de exemplu în cartea lui Hoerkins și Hails ( ) TRANSFER ȘI ROTIRE ÎN SPAȚIU TRIDIMENSIONAL Dacă fiecare punct P(x, y, z) este mapat pe un punct P'(x', Y'" z) în conformitate cu ecuaţiile x' = x + a l + arctan (v /vp dacă ѵ^O l/ dacă Vj = și v > Zl/ dacă Vj = și v COS Ѳ, V =p sin sin Ѳ, V =p cos ) - Acum să rotim axele x' și z' în jurul axei y' cu un unghi R Despre -șina O DESPRE cosa Rotația reală în jurul vectorului ѵ prin unghiul a poate fi acum efectuată ca o rotație în jurul axei z" Din ecuația ( ) obținem y'" z"']= O" y" z"]Rv cosa șina Лѵ= -șina cosa aproximativ În acest moment, relația [x'" y"' z"']= [x y z]R~lR'lRv Orez Coordonate sferice TRANSFER ȘI ROTIRE ÎN SPAȚIU TRIDIMENSIONAL Din păcate, coordonatele x'", y"', z'" se referă la cel mai recent sistem de coordonate, în timp ce ele trebuie exprimate în sistemul original Să notăm aceste coordonate în sistemul original cu x*, y*, z* sistem de matrici inversate R~x si R~x (care vor coincide cu matricele R si R) in ordine inversa pentru a transforma punctul x'", y'", z Y z : [x* y* z*] = [x'" y'" z"'] zbor/îz Aceasta înseamnă că o rotație completă în jurul vectorului ѵ după unghiul a este calculată prin următoarea formulă: [x* y* z*]=[x"' y'" Unde ' cos Ѳ -sin Ѳ sin# cos# V cos cos^" cos# sin# O' R= -sin# cos# O [o Pentru utilizare ulterioară scriem " ■'V'YALL-I' r r r r r r J ( , ) Până acum, am discutat soluția problemei rotației cu privire la vector atașat la punctul de origine al sistemului de coordonate O Acum trebuie să eliminăm această ultimă limitare și să postăm rezolvați problema determinării rotației în jurul vectorului, al cărui început este situat în orice punct arbitrar A(aț, a , л ) Capitolul INSTRUMENT GEOMETRIC Pentru a face acest lucru, vom folosi vectorul ѵ pentru a calcula matricea R în ecuația ( ) în același mod ca înainte Apoi, trebuie să urmați următorii trei pași: Revenind la ecuația ( ), vom efectua transferul de la un punct dat la punctul de origine O, folosind un coordonatele și următoarea matrice: O' -a, -a -a Acum ne putem roti în jurul axei care trece prin O, ca înainte, dar matricea R din ecuația ( ) trebuie extinsă într-un mod trivial, astfel încât să poată fi utilizate coordonate omogene GP G G G G G G G G OOO ooo Aplicați transformarea inversă a pasului folosind matricea ' O' a a a - După aceea, matricea de rotație generalizată este calculată ca rge# '- *T si poate fi folosit astfel: [x* y* z* ] = [x y z ] tfGEN EXERCIȚII Explorați efectul alegerii a = / = y " în ecuația ( ) din secțiunea și scrieți un program pentru a mapa o serie de puncte dintr-un cadran la un triunghi dat Obțineți imagini triunghiulare și punctuale ( , ), ( , ), ( , ), ; ( , ), ( , ), ( , ), și așa mai departe sub formă grafică EXERCIȚII Scrieți un program pentru o rotație generalizată Datele de intrare pentru program vor fi: ) coordonatele a două puncte care definesc vectorul ѵ = AB cu originea în punctul A; ) unghiul a Rotirea se va efectua în jurul vectorului ѵ prin unghiul a; ) o serie de puncte de rotit Rezultatul programului ar trebui să fie coordonatele punctelor rotite Fiecare punct este descris de un trio de coordonate într-un sistem ortogonal Deoarece mijloacele de ieșire grafică a obiectelor tridimensionale nu au fost încă luate în considerare, rezultatele programului trebuie să fie afișate în formă numerică Găsiți un poligon care nu poate fi procesat corect de programul POLY TRIA și dezvoltați un algoritm mai general Sugestie: Dacă un triunghi taie o parte dintr-un poligon, atunci niciunul dintre vârfurile acestui poligon nu ar trebui să cadă în interiorul triunghiului capitolul IMAGINI DE PERSPECTIVĂ INTRODUCERE Pe fig Figura prezintă o imagine bidimensională a unui cub împreună cu câteva linii suplimentare În această imagine, segmentele de linie AB și AD nu sunt paralele cu marginea de jos sau de sus a hârtiei, așa că s-ar putea argumenta că nu sunt orizontale Cu toate acestea, ele denotă muchiile orizontale ale cubului ABCDEFGH în spațiul tridimensional, deci în principiu pot fi numite orizontale Din același motiv, se poate argumenta că cele două segmente AB și DC sunt paralele, presupunând implicit că se află în spațiu tridimensional În această terminologie, liniile orizontale paralele se întâlnesc în așa-numitul punct de fugă Toate punctele de fugă se află pe o singură linie dreaptă, care se numește linia orizontului Rețineți că linia orizontului și punctul de fuga sunt caracteristici ale imaginii și nu există cu adevărat în spațiul D Timp de secole, acest concept a fost folosit de artiști pentru a crea imagini realiste ale obiectelor D Astfel de imagini sunt de obicei numite perspectivă Invenția fotografiei a oferit un mod nou (și mai ușor) de a forma imagini în perspectivă Existent- Orez Puncte de fuga la orizont INTRODUCERE Există o analogie strictă între aparatul foto folosit în fotografie și ochiul uman Ochiul nostru este un instrument foarte complex, iar camera este cea mai simplă imitație a sa Dar în discuția care urmează, cuvântul ochi poate fi bine înlocuit cu cuvântul cameră, dacă vrem să subliniem că este de dorit să obținem o copie bidimensională pe hârtie Evident, imaginea va depinde de poziția ochiului De o importanță deosebită este distanța dintre ochi și obiect, deoarece "efectul de perspectivă" va fi invers proporțional cu această distanță Dacă ochiul este foarte aproape de obiect, atunci obținem un efect de perspectivă puternic, ca în Fig (a) Aici puteți vedea clar că continuarea imaginilor liniilor paralele din imagine se intersectează Pe de altă parte, dacă ochiul este departe de obiect (comparativ cu dimensiunea obiectului), atunci liniile paralele ale obiectului vor apărea paralele și în imagine Acest lucru este prezentat în fig ( ) Orez Locația ochiului: (a) aproape de obiect; (b) - departe de obiect Pe lângă metodele clasice și fotografice, există o metodă de obținere a imaginilor în perspectivă bazată pe geometria analitică Cititorul este deja familiarizat cu reprezentările punctelor din spații bidimensionale și tridimensionale prin coordonatele lor (x, y) și respectiv (x, y, z) Dacă este necesar să se obțină o proiecție în perspectivă, se precizează un număr mare de puncte P(x, y, z) aparținând obiectului, pentru care coordonatele punctelor imaginii P'(X, Y) din imagine sunt să fi calculat Pentru a face acest lucru, trebuie doar să convertiți coordonatele punctului P din așa-numitele coordonate mondiale (x, y, z) în coordonatele ecranului (X, Y) ale proiecției sale centrale?' Vom presupune că ecranul este situat între obiect și ochiul E Pentru fiecare punct P al obiectului, o linie dreaptă PE intersectează ecranul în punctul P' Capitolul IMAGINI DE PERSPECTIVĂ Această mapare este efectuată în mod convenabil în doi pași Prima etapă va fi numită o transformare a vederii - punctul P rămâne la locul său, dar sistemul de coordonate mondial intră în sistemul de coordonate a vederii A doua etapă se numește transformarea perspectivei Aceasta este transformarea exactă a punctului P în punctul P', combinată cu tranziția de la sistemul de coordonate a vederii D la sistemul de coordonate a ecranului D: Coordonatele lumii (xw, yw, zw) Vizualizați transformarea Vedeți coordonatele (he, ye, ze) Transformarea perspectivei j Coordonatele ecranului (X, Y) VEZI CONVERSIUNEA Pentru a efectua transformări de vedere, trebuie specificat un punct de observație care coincide cu ochiul și un obiect Este de dorit ca sistemul de coordonate mondial să fie corect Va fi convenabil dacă originea coordonatelor sale este situată undeva în apropierea centrului obiectului, deoarece obiectul este observat în direcția de la E la O Să presupunem că această condiție este îndeplinită În practică, aceasta înseamnă că poate fi necesară o anumită transformare a coordonatelor, care constă în scăderea poziției punctului central al obiectului din valorile originale ale coordonatelor Această transformare foarte simplă va fi inclusă în program, dar nu va fi scrisă în formă matematică Fie dat punctul de observație E în coordonate sferice /), Ѳ, -l Din ecuația ( ), obținem matricea: ' Yax \u d costyj-l) sinty\u e -l) [o -Z P (^-L) costyj-l) :G = -cos -sin# sin Ѳ -cos Ѳ ( , ) Noile axe sunt prezentate în fig Schimbarea direcției axei x Pe fig axele y și z au orientarea corectă, iar axa x trebuie să fie îndreptată în direcția opusă Prin urmare, este necesară o matrice pentru a efectua transformarea x' = -x, adică myz= - sau o ( , ) După această transformare finală, obținem sistemul de coordonate a vederii prezentat deja în Fig VEZI CONVERSIUNEA Calculăm matricea de mapare V ca produs de matrice V= TR*R*My* ( , ) unde R* este folosit pentru o matrice * obținută prin extinderea unei matrice * R prin adăugarea unui al patrulea rând și a unei a patra coloane care conține numerele , , , în această ordine Produsul matriceal nu este comutativ (adică în cazul general AB * B A), ci este asociativ, deci ecuația ( ) poate fi rescrisă ca i-t(/glmu/ Deci puteți lucra cu matrice x cât mai mult posibil Elaborarea ulterioară a problemei înmulțirii matricelor este asociată cu utilizarea notației cos F \u d a cos Ѳ \u d s / Och , n, ( , ) sin -cos Z Mai trebuie să determinăm distanța dintre punctul de observare E și ecran În linii mari, avem relația: dimensiunea imaginii dimensiunea obiectului d - p care rezultă din asemănarea triunghiurilor EP'Q' şi EPQ din fig Orez Dimensiunile imaginii și obiectelor Capitolul IMAGINI DE PERSPECTIVA De aici ajungem dimensiunea imaginii ( ) n dimensiunea obiectului Această expresie este valabilă în egală măsură pentru dimensiunile orizontale și verticale Ar trebui interpretat mai mult ca un mijloc de a estima o valoare adecvată decât o prescripție precisă, deoarece un obiect tridimensional poate avea o formă foarte complexă și nu este întotdeauna clar ce dimensiuni ar trebui incluse în această ecuație Puteți introduce o estimare aproximativă a dimensiunilor unui obiect pe baza maximului de lungime, lățime și înălțime Pe baza evaluării ecuației ( ), putem concluziona că dimensiunea imaginii ar trebui să fie ceva mai mică decât dimensiunea ecranului Modalități mai sofisticate de estimare a dimensiunilor dorite ale imaginii vor fi descrise în Secțiunea Restul acestei secțiuni este dedicat unei analize mai detaliate a conceptului general, în special în ceea ce privește punctele de fugă și orizont Pe fig arată punctele de vedere E și ecranul ABCD Ca și până acum, privirea este îndreptată din punctul E către punctul Q Liniile drepte AF, BG, EH sunt paralele și le vom considera orizontale Să presupunem că există un plan care trece prin liniile paralele EH și BG Acest plan intersectează planul ecranului într-o linie dreaptă HV Astfel, fiecare punct P de pe linia dreaptă BG va avea proiecția sa centrală?" situată pe dreapta BH, cu condiția ca punctul E să fie centrul proiecției Fie ca punctul P să se îndepărteze de punctul B la infinit, atunci proiecția sa P' va fi Orez Punctul de fuga H TRANSFORMĂRI DE PERSPECTIVĂ se apropie de punctul H Aceasta înseamnă că H este punctul de fugă al dreptei care trece prin punctele B și G În ceea ce privește geometria proiectivă, punctul H este proiecția unui punct la infinit situat pe dreapta BG După cum sa menționat în Secțiunea , liniile paralele se intersectează într-un punct la infinit, deci punctul H va fi, de asemenea, o proiecție a unui punct la infinit situat pe o linie dreaptă AF Dacă luăm o linie dreaptă cu o direcție diferită, dar și una orizontală, de exemplu, o linie dreaptă BF, atunci astfel de Orez Punctul de fuga F pentru linii verticale Capitolul IMAGINI DE PERSPECTIVA linia va avea, de asemenea, un punct de fuga situat pe linia dreaptă CD Se găsește ca punct de intersecție al dreptei CD și al dreptei care trece prin punctul E și paralel cu această dreaptă orizontală aleasă Linia dreaptă CD este linia orizontului Fiecare punct! pe linia orizontului este punctul de fugă al tuturor dreptelor paralele cu dreapta EJ Punctul de fuga nu este doar pentru liniile orizontale Pe fig sunt trasate linii drepte verticale CA și DB Punctul de fuga pentru ei va fi F Acesta este punctul în care linia dreaptă verticală care trece prin punctul E traversează ecranul Segmentele de linie dreaptă CA și DB au proiecțiile CA' și DB' care nu sunt paralele Un efect de perspectivă strict pentru liniile drepte verticale nu este întotdeauna acceptabil; imaginile în care liniile verticale apar aproape verticale sunt adesea preferate Acest lucru se datorează faptului că suntem mai obișnuiți cu o linie de vedere orizontală sau aproape orizontală Unii chiar se simt amețiți dacă privesc departe în jos! Artiștii folosesc "pseudo-perspectivă" în picturile lor, unde liniile verticale sunt desenate exact verticale, chiar dacă direcția de vedere nu este orizontală În acest caz, imaginea rezultată diferă de cea cu adevărat vizibilă, dar, în ciuda unei oarecare curiozități, pare destul de plauzibilă Un exemplu este Fig ( ) În Secțiunea , vom reveni la acest fenomen și vom arăta că programul nostru poate genera o astfel de pseudo-promițătoare Orez (a) - perspectiva; (b) - pseudoperspectivă TRANSFORMĂRI DE PERSPECTIVĂ ABCDEFGH J Orez Zece zaruri paralele cu ecranul imagini (vezi Figura (a)) În general, însă, se recomandă alegerea unui punct de observare nu foarte aproape de obiect, mai ales în cazurile în care unghiul = , c/= și fig ( ) cu valori p = , * , ) dacă (cod) dw(x-xO, y-yO, z-zO); else mv(x-xO, y-yO, z-zO); endgra(); } skipf(fp) FIȘIER *fp; { while (getc(fp)!- '\n'); } coeff(rho, theta, phl) float rho, theta, phi; { float th, ph, costh, slnth, cosph, sinph, factor; factoratan( , J/ , ; /* Unghiuri în radlani: */ /* Unghiuri în radiani */ factor-th-theta*; ph-phi*factor; cost-cos(th); slnth-sin(th); cosph-cos(ph); sinph-sin(ph); /♦ Elemente ale matrlx V, vezi Eq ( , ); */ /* Elemente Matrix V, vezi ( ): */ v -~sintet; v -cosph*cost; v -sinph*cost; DESENARE MODELE DE SÂRMĂ v -cost; v -cosph*slnth; v -slnph*slnth; v -slnph; v -cosph; v -rho; } mv(x, y, z) float x, y, z; { floatXY; perspectlve(x, y, z, &X, &Y); mutare(X, Y); } dw(xyz) float x y, z; { floatXY; perspectlve(x, y, z &X, &Y); trage (X, Y); } perspectlve(x, y z, pX pY) float x, y, z *pX, *pY; { float x ye, ze; /* Coordonatele ochiului, calculate ca în Ec ( ): */ /* Vedeți coordonatele calculate de nr ( ): */ xe-v *x + v l*y; ye - v *x + v *y + v *z; ze - v *x + v *y + v *z + v ; /* Coordonatele ecranului calculată ca InEqs ( ) și ( ): */ /★ Coordonatele ecranului calculate din ( ) și ( ) */ *pX - screen dlst*xe/ze + c ; *pY - screen dlst*ye/ze + с ; } Funcția fscanf ia o valoare nepozitivă atunci când o încercare de citire a datelor eșuează, adică atunci când se ajunge la sfârșitul fișierului Prin urmare, bucla de citire continuă atâta timp cât funcția fscanf returnează o valoare pozitivă Utilizatorul programului GPERS trebuie să seteze distanța până la ecran Din mai multe motive, această soluție este nesatisfăcătoare: Mărimea imaginii este importantă pentru utilizator, nu distanța până la ecran Parametrul și dimensiunea "obiectului" din formula ( ) are o semnificație foarte vagă, deci este dificil de calculat valoarea exactă pentru parametrul d Uneori este de dorit să definiți o zonă de ieșire dreptunghiulară, care ocupă doar o parte a ecranului, în care imaginea trebuie introdusă în întregime În acest caz ar fi foarte Capitolul IMAGINI DE PERSPECTIVA este incomod să convertiți dimensiunile zonei de ieșire la distanța până la ecran Există mai multe modalități de îmbunătățire a programelor în aceste domenii, și anume: Tăierea unui obiect tridimensional de-a lungul unei piramide al cărei vârf coincide cu punctul de observație E, iar baza este o anumită fereastră specificată în coordonatele lumii O descriere a unei astfel de metode, care nu va fi folosită aici, poate fi găsită în Newman și Sproul ( ) Tăierea imaginii în spațiu bidimensional în funcție de zona de ieșire specificată Nu vom folosi nici această metodă Selectarea automată a dimensiunii și poziției astfel încât întreaga imagine să fie inclusă în zona de ieșire specificată Am folosit deja această metodă pentru obiecte bidimensionale în Secțiunea , unde programul GENPLOT a acționat ca un post-procesor Același program poate fi folosit în acest caz Mai mult, vom putea folosi ultimele treisprezece linii ale programului CURVGEN din Secțiunea pentru a scoate datele de trasare în fișierul A SCRATCH Vom folosi a treia metodă, folosind două fișiere și două programe, așa cum se arată în Fig În programul GPERSF, valoarea lui d nu este introdusă, ci pur și simplu setată la d = Această valoare practic nu este necesară, deoarece dimensiunile imaginii calculate de programul GENPLOT depind doar de dimensiunile zonei de ieșire specificate Amintiți-vă că programul GENPLOT cere utilizatorului limitele zonei de ieșire În această etapă, cititorul se poate întreba dacă este posibil să scape de necesitatea precizării punctului obiect O, deoarece punctul O a fost descris ca punct central al obiectului și se pare că poziția unui astfel de punct poate fi ușor de calculat automat Cu toate acestea, există două motive pentru care ar trebui să lăsați punctul O ca punct definit de utilizator Aceste motive vor fi discutate mai jos în Secțiunea Următorul program GPERSF preia același fișier de intrare ca și programul GPERS, dar distanța ecranului de la sfârșitul celei de-a doua linii de intrare este ignorată și poate fi omisă DESENARE MODELE DE SÂRMĂ ѵ Ieșire grafică Orez Diagramă pentru programele GPERSF și GENPLOT I* GPERSF: Program general pentru PERSpectlve, care produce */ /* Fișierul de ieșire A SCRATCH pentru a fi citit de GENPLOT */ /* Program generic de imagine în perspectivă */ /* generează un fișier de ieșire A SCRATCH pentru a fi */ /* citit de GENPLOT */ tfinclude "math h" tflnclude "stdio h" plutitor v , v , v , v , v , v , v , v , v ; maln(argc, argv) Int argc; char *argv[ ]; { float rho, theta, phi, xO, yO, zO, x, y, z; FIȘIER *fpln, *fpout; struct { float X, Y; intcode; }s; If (argc - ) { printf /* "Fără fișier de intrare glven" */ ("Nici un fișier de intrare specificat"); ieșire( ); } if (fpin-fopen(argv[ ],"r"), fpin- -NULL) { printf /* "Fișierul nu există" */ ("Fișierul %s nu există", argvflj; exlt( ); } ON Capitolul IMAGINI DE PERSPECTIVA fscanf(fpln, "%f %f %r, &x , &y , &zO); sklpf(fpin); fscanf(fpln, "%f %f %f", &rho, &theta, &phi); coeff(rho, theta, phl); fpout-fopen ("a scratch", "wb"); în timp ce (sklpf(fpln), fscanf(fpin, "%f %f %f %d" &x, &y, &z, &s code)> ) { perspectiva(x-xO, y-yO, z-zO, &s X, &s Y); fwrlte(&s, sizeof s, , fpout); } fclose(fpout); } sklpf(fpln) FIȘIER *fpln; { while (getc(fpln)!- '\n'); } coeff(rho, theta, phl) float rho, theta, phl; { float th, ph, costh, slnth, cosph, slnph, factor; factor-atan( , )/ , ; /* Unghiuri în radlani: */ /* Unghiuri în radiani */ th-theta*factor, costh-cos(th); cosph-cos(ph); ph-phi*factor; slnth-sln(th); slnph-sln(ph); /* Elementele matrlx V, vezi Eq ( ): */ /* Elementele matricei V, vezi ( ): */ v -slnth; v -cost; v -cosph*cost; v -slnph*cost; v -cosph*slnth; v -slnph*slnth; v -slnph; v -cosph; v -rho; } perspect(x, y, z, pX, pY) float x, y, z, *pX, *pY; { float xe, ye, ze; /* Coordonatele ochiului, calculate ca în Ec ( ) */ /* Vedeți coordonatele calculate de nr ( ) */ xe- v *x + v *y; ye - v *x + v *y + v *z; ze - v *x + v *y + v *z + v ; /* Coordonatele ecranului, pot fi ajustate de GENPLOT */ /* Coordonatele ecranului recalculate în programul GENPLOT */ *рХ - xe/ze; *pY - ye/ze; } Funcția de perspectivă din acest program a fost derivată din funcția de perspectivă din programele anterioare CUBE și GPERS prin înlocuirea ecranului jiists și cl s c e (Alte valori vor schimba conținutul fișierului A SCRATCH, dar imaginea generat de programul GENPLOT va fi exact același!) DIRECȚIA DE OBSERVARE DIRECȚIA OBSERVAȚIEI, INFINITATE, LINII VERTICALE Să folosim programul GPERSF pentru a discuta câteva aspecte noi interesante Să presupunem că obiectul este foarte lung în direcția x, cum ar fi fasciculul din Fig , ale căror dimensiuni sunt x x Un aspect interesant al acestui exemplu este precizarea punctului O, pe care trebuie să-l precizăm pentru a determina direcția de observare a EO Până acum, punctul O a fost ales în centrul obiectului Cu toate acestea, de fapt, este necesar un astfel de punct O, a cărui imagine va fi situată în centrul imaginii Adesea, acest lucru nu este semnificativ, dar uneori afectează foarte mult imaginea De exemplu, în fig Punctul O' se dovedește a fi în mijlocul segmentului Q'U', deși punctul O corespunzător acestuia nu este ales în mijlocul segmentului original QU Dacă ar fi dat punctul de mijloc al segmentului QU, atunci direcția de observație s-ar dovedi a fi complet diferită Să revenim la fasciculul prezentat în fig Evident, punctul O nu trebuie ales în mijlocul fasciculului, ci mult mai aproape de ochi Acest lucru este posibil deoarece în programul nostru punctul obiect este dat de utilizator și nu calculat ca centru al obiectului Liniile reale fără sfârșit atunci când înfățișați un peisaj în imagine sunt obținute sub formă de segmente Capitolul IMAGINI DE PERSPECTIVA linii drepte de lungime finită Pentru obiecte infinite, nu are rost să vorbim despre centru, deși direcția de observare a EO există Prin urmare, este necesar să se determine poziția punctului "central" O al obiectului Grinda din fig , desigur, nu este un peisaj în sensul exact, dar lungimea lui sugerează existența infinitului; dacă acest lucru nu este complet clar, atunci pentru a determina lungimea sa, puteți specifica un număr mult mai mare decât Pentru acest fascicul, selectam punctul O (- , , ) Cu excepția literelor P, Q, , toată această imagine a fost desenată folosind programele GPERSF și GENPLOT, intersecția fasciculului cu planul x = r este indicată prin literele ABC, iar intersecția cu x = - plan cu literele A'B'C ' Prin urmare, punctul O se află în planul ABC, iar planul A'B'C' este situat în centrul fasciculului Setul complet de date de intrare pentru obținerea fig definite de lista: - (obiect punctul O) (ro, tega, fi) (deplasare la punctul Q) (desenarea unui segment QR) (desenarea unui segment RS) (desenați segmentul de linie SP) (segment de desen PQ) - (segment de desen QU) - (tras linie UV) - (desenați segmentul de linie VW) DIRECȚIA DE OBSERVARE (segment de desen WS) ' (deplasare la punctul R) - - (desenați segmentul de linie RV) - (deplasare la punctul A) - (desen segment AB) - (segment de desen BC) - (deplasare la punctul A') - (segment de desen A'B') - (desen segment B'C') Un alt aspect care necesită o atenție specială este aceasta este reprezentarea liniilor verticale, menționată pe scurt în secțiunea Pe fig În Figura , am văzut că liniile verticale proiectate se întâlnesc în punctul de fugă F De asemenea, am comparat reprezentările cuburilor din Figura Dacă punctul O este situat în centrul cubului, atunci programul nostru nu va putea forma imaginea din Fig ( ) Dar punctul O poate fi plasat peste cub astfel încât direcția de observație EO să fie orizontală Apoi ecranul va ocupa o poziție verticală, deoarece este perpendicular pe direcția de observare Aceasta înseamnă că marginile verticale ale cubului vor fi paralele cu ecranul Prin urmare, proiecțiile lor în imagine vor fi, de asemenea, exact verticale Lucrul curios aici este că, deși linia de observație este orizontală, vom vedea fața superioară a cubului ca de sus Poza rezultată pare complet naturală, atâta timp cât nu exagerăm și nu alegem punctul de vedere prea sus Pe fig arată ambele astfel de cazuri Liniile din fig (a) au fost reprezentate grafic de programele GPERSFh GENPLOT pe baza următoarei liste de intrări: (punctul O) (ro, tega, fi) - (deplasarea în punctul P) - - (segment de desen PQ) - (segment de desen QU) (segment de desen UT) - (trageți segmentul TR) - - (desenați segmentul de linie PR) - (segment de desen RW) Capitolul IMAGINI DE PERSPECTIVA (desenați segmentul de linie WT) - (deplasare la punctul W) - - (segment de desen WV) - (segmentul VU desenat) (axa x) (axa y) (axa z) Orez ( ) se obține prin alegerea unui punct O mult mai mare și anume, cu coordonatele ( , , ) în loc de ( , , ) Aceasta este o imagine complet nerealistă a cubului În același timp, Fig (a) este perfect acceptabil Unii chiar îl preferă în detrimentul cuburilor de orez (a) și fig (a), unde liniile verticale au un punct de fugă Dacă cubul este desenat manual în perspectivă, atunci de obicei liniile verticale rămân verticale, ca în fig (a) Alte reprezentări sunt mult mai dificil de realizat manual Cu ajutorul unui computer, acestea sunt la fel de ușor de executat și utilizatorul le poate alege pe cele care îi plac Orez (a) punctul obiectiv este ceva mai înalt decât cubul; (b) - punctul obiectului este foarte sus deasupra cubului (b) EXERCIȚII cometariu Programele din acest capitol au scopul de a explica mai degrabă principii decât aplicații practice Capitolele și vor discuta mai multe programe practice EXERCIȚII În următoarele exerciții, marginile invizibile nu trebuie desenate, pentru aceasta puteți folosi metoda descrisă la sfârșitul paragrafului Scrieți un program pentru a desena pași ca cei din Fig Numărul de pași (l) trebuie să fie variabil (aici n în ) Scrieți un program pentru a desena o piramidă de cuburi ca cea prezentată în Fig Orez Piramidă Capitolul IMAGINI DE PERSPECTIVA Numărul n, care determină înălțimea piramidei în unități de cuburi, trebuie citit de program (aici n - ) Scrieți un program pentru a desena un obiect ca cel prezentat în fig Acesta este un cub din care au fost tăiate mai multe cuburi mai mici; dimensiunile lor sunt /l-th din cubul original Numărul n trebuie să fie variabil (aici l = ) capitolul ELIMINAREA LINIILOR ASCUNSE ANALIZA PERFORMANȚEI Programele simple nu sunt întotdeauna mai bune decât cele complexe În unele probleme, algoritmii simpli sunt prea lenți și necesită algoritmi mai puternici, uneori implementați ca programe extrem de complexe Un exemplu ar fi sortarea Programele descrise în capitolul pentru generarea desenelor în perspectivă au fost relativ simple și rapide Evident, complexitatea unor astfel de programe va crește semnificativ dacă va fi necesară eliminarea automată a părților invizibile ale segmentelor Odată cu complexitatea sarcinii în sine, apar inevitabil și alte probleme atunci când programul trebuie să ruleze cât mai repede posibil Pe lângă asigurarea eficienței, este necesar să se satisfacă cerințele de generalitate, toleranță la erori și ușurință în utilizare, ceea ce duce în cele din urmă la potențiala complexitate a programului S-a acordat atât de multă atenție acestor patru aspecte ale calității încât cititorul ar putea crede că acum vom lua în considerare doar programe foarte complexe De fapt, este mai atractiv să începi cu un program oarecum simplificat decât cu unul rapid, dar de neînțeles Prin urmare, începem cu un program cu complexitate de ordin n Aceasta înseamnă că timpul de calcul va fi aproximativ proporțional cu n , unde n este numărul de segmente de linie dreaptă sau poligoane de desenat Programul va funcționa suficient de rapid atunci când desenați obiecte simple care conțin câteva sute de fețe Datele care descriu obiectul vor fi stocate în matrice de dimensiuni fixe De asemenea, va limita complexitatea obiectelor Pe lângă aceste restricții, restul programului va fi destul de general: în principiu, oricare Capitolul ELIMINAREA LINIILOR ASCUNSE Foste obiecte finite care au doar un număr finit de fețe plane În ceea ce privește toleranța la erori și confortul utilizatorului, nu vom cere ca toate erorile, atunci când apar, să fie însoțite de mesaje clare și ca obiectele să fie definite în cel mai convenabil mod Ulterior acest program va fi îmbunătățit, mai ales din punct de vedere al eficienței Problema de eliminare a liniilor ascunse este implementată de algoritmi extrem de intensivi din punct de vedere computațional și a fost întotdeauna o problemă bună pentru matematicienii programatori INTRARE ȘI REPREZENTARE INTERNĂ Marginea unui obiect poate fi acoperită complet sau parțial de una sau mai multe fețe ale aceluiași obiect Fiecare muchie este segmentul de capăt al unei linii drepte Un obiect poate consta din mai multe părți, nu neapărat interconectate Pe fig Figura prezintă un tetraedru și un cub cu vârfuri etichetate , , , în loc de literele A, B, C, Punctele etichetate , , vă permit să desenați semi-axe de coordonate pozitive Orez Tetraedru și cub INTRARE ȘI REPREZENTARE INTERNĂ Să scriem un program care va citi datele care descriu punctul de observație și obiectele Pentru orice punct de vedere (real), rezultatul grafic va fi o imagine a obiectului Spre deosebire de capitolul , vom presupune că obiectul este opac Ca și în Secțiunea , utilizatorul trebuie să furnizeze două linii de date de intrare specificând punctul O în sistemul de coordonate carteziene și coordonatele sferice p, pentru punctul de vedere Pentru fig vor fi următoarele rânduri: (centrul obiectului O) (rho, theta, phi) Nu vom acorda prea multă atenție confortului utilizatorului aici, ci mai degrabă introducem conceptul de codificare a vârfurilor: deoarece fiecare vârf va fi necesar de mai multe ori, ar fi foarte nepractic să specificați coordonatele sale de mai multe ori În acest exemplu, pentru fiecare vârf îi vom seta numărul de la la , împreună cu valorile coordonatelor de-a lungul axelor x, y, z Sfârșitul acestei părți a informațiilor de intrare va fi indicat printr-o linie separată cu simbolul # în prima poziție: (x = , y = , z = ) (xi = , yi = , zl = ) (etc ) # Ca și în Secțiunea , punctul O al obiectului este originea sistemului de coordonate mondial utilizat în program (rețineți diferența dintre litera O și numărul ) Deoarece are coordonate ( , , ), atunci în cadrul programului coordonatele punctelor de la la vor fi reduse cu aceste valori Deci, de exemplu, punctului i se vor atribui valori de coordonate (- , , , ) în locul celor date ( , , ) Fiecare față a unui obiect poate fi descrisă ca orice zonă finită a planului cu limite sub formă de segmente de linie dreaptă Un exemplu de astfel de regiune este poli Capitolul ELIMINAREA LINIILOR ASCUNSE gon, în care sunt permise și găuri Pentru a simplifica programul, utilizatorul trebuie să împartă el însuși astfel de zone în triunghiuri Aceste triunghiuri vor fi apoi folosite pentru a determina dacă acopera segmente de linie dreaptă Din motive care vor fi explicate mai jos, numerele vârfurilor din fiecare triunghi sunt enumerate în ordinea inversă acelor de ceasornic atunci când sunt privite din exteriorul obiectului Fiecare față a cubului este împărțită în două triunghiuri Această parte a intrării este din nou terminată cu simbolul # (pe o linie nouă) (fețele triunghiulare ale unui tetraedru) (fața frontală a cubului) (partea dreaptă) (linia superioară) (fața din spate) (partea stângă) (jos) # Acum trebuie să definim fiecare margine a obiectului Deși aceste margini sunt deja cunoscute ca laturile triunghiurilor, ele ar trebui descrise din nou În primul rând, nu toate laturile triunghiurilor sunt margini ale obiectului și, în al doilea rând, este de dorit să se poată desena segmente suplimentare care nu aparțin unui corp solid Pentru ilustrare, desenăm părți din coordonatele pozitive ale semiaxelor, așa cum se arată în Fig Destul de curios, în acest exemplu, numărul liniilor de intrare nu va crește, deoarece marginile obiectului care se află pe axele de coordonate nu trebuie să fie redefinite Prin urmare, avem linii de intrare: INTRARE ȘI REPREZENTARE INTERNĂ (axa x) (axa y) (AXA Z) (tetraedru) (cub) Programul va citi toate intrările din fișierul al cărui nume este transmis ca argumente de program (argc și argv), așa cum este descris în Secțiunea Coordonatele fiecărui vârf sunt stocate în tabloul VERTEX, ale cărui elemente sunt structuri care conțin trei câmpuri x, y, z Coordonatele utilizatorului specificate sunt mai întâi convertite în coordonate de sistem, a căror origine este situată în punctul obiect O Coordonatele lumii interioare, la rândul lor, sunt convertite în coordonatele vizualizării folosind transformări de vizualizare (vezi Secțiunea ) Și aceste coordonate de vizualizare sunt scrise în tabloul VERTEX* i O VERTEX[i] X y ZII După cum știm deja, punctul E este originea sistemului de coordonate a vederii, iar direcția de observare coincide cu direcția axei z Prin urmare, toate valorile VERTEX [i ] z trebuie să fie pozitive Această condiție este verificată în program, iar execuția programului se oprește dacă nu este îndeplinită Lista de triunghiuri este stocată în tabloul TRIANGLE* j TRIANGUL C] I A B C a b c h eu voi DESPRE Capitolul ELIMINAREA LINIILOR ASCUNSE Împreună cu numărul de vârfuri pentru fiecare triunghi se scriu coeficienții a, , c, h ai ecuației planului ax + by + cz = h ( ) în care se află triunghiul Desigur, ar putea fi calculate de fiecare dată când sunt necesare, dar ar fi o pierdere de timp, deoarece sunt solicitate destul de des (Aici al patrulea coeficient este notat cu h în loc de c/, deoarece această literă este deja folosită pentru a indica distanța până la ecran ) Coeficienții a, b, c sunt aleși astfel încât condiția să fie îndeplinită a + b + c = și A > O Atunci ecuația ( ) poate fi scrisă ca produs scalar n • x = A Unde P \u d [ ) { r = sqrt(a* a + b* b + c* c); a = a!r\ bb!r\ cc/r\h^h/r\ } altfel { /* Triunghiul ABC poate fi ignorat din următoarele motive */ } Un singur calcul al coeficienților a, , c, A, în loc să-i determine la fiecare verificare a vizibilității segmentului, poate reduce semnificativ timpul de calcul Dacă programul chiar trebuia să fie cât mai simplu posibil, atunci ar fi necesar să memorezi toate triunghiurile pe măsură ce au fost introduse Cu toate acestea, acele triunghiuri care se află în spate nu trebuie să fie salvate și pot fi ignorate Se consideră triunghiul din fig În imagine, este închis de triunghiul Deoarece ultimul triunghi acoperă o parte a segmentului , putem ignora faptul că triunghiul anterior face același lucru Fețele din spate sunt acoperite de fețele vizibile Deși fețele din spate pot acoperi unele puncte din ochi, aceste puncte sunt, de asemenea, acoperite de margini vizibile De aceea, fețele din spate pot fi ignorate Cel mai simplu mod de a identifica fețele din spate se bazează pe orientarea vârfurilor Dacă ne uităm la fața din Fig în spațiul D din exterior (adică din partea semiaxei negative x), atunci ordinea de parcurgere a vârfurilor la intrare va fi în sens invers acelor de ceasornic, deoarece o astfel de orientare era necesară la definirea Capitolul ELIMINAREA LINIILOR ASCUNSE nii a secvenței de intrare Cu toate acestea, în fig ordinea bypass-ului corespunde mișcării acelui orelor Aceasta înseamnă că în imagine această față este vizibilă prin corpul obiectului, și nu din exterior Astfel, fata a revenit Aplicarea acestei metode de determinare a poziției fețelor pentru alte triunghiuri din fig poate servi ca un exercițiu bun Referindu-ne la sfârșitul secțiunii , se constată că în acest caz este necesar să se găsească valoarea determinantului D= hk *v xs UA unde ХА, Уд, Хв, Ув, Хс, Ус sunt coordonatele de ecran ale vârfurilor А, В, С ale triunghiului Triunghiul ABC va fi înapoi dacă D = XA Wx ZA XB ZB xs Ys zc ! {zK ZB zcj = b/(zA zB zc) ( , ) Deoarece valorile lui zA, zfî, zc sunt întotdeauna pozitive, atunci D are același semn ca și valoarea lui A calculată în programul anterior Dacă A \u d , atunci planul ABC trece prin punctul de origine E, care coincide cu punctul de observație În acest caz, triunghiul ABC nu va acoperi niciun segment și nu trebuie să fie stocat în matricea TRIANGLE Când A sau, cu alte cuvinte, dacă D > Cu ajutorul produsului încrucișat, se poate asigura că triunghiul ABC este înapoi, din condiția A z-zA+/z(zB-zA) Pentru valorile lui /i de la la , punctul se află între A și B Ecuația planului EPQ X y *r Ur XQ Să rescriem această ecuație sub forma KuX + K y + K^z = z ZP ZQ = Unde ^ = XQZP XPZQ Wq'M" Combinând ecuațiile ( ) și ( ), obținem KlxK + K yA + K z\ " Xp) + " Y+ K (ZB " Un segment de linie dreaptă PQ și planul piramidei EAB au un punct comun dacă și numai dacă ( , ) ( , ) ( , ) epsl x a x> a + epsl x>°a x>° a- epsl x = = fabs(x) epsilon x x > epsilon x>= x >= -epsilon Rețineți că a doua parte a acestui tabel poate fi derivată din prima Aceste teste vor fi utilizate în programul HIDLIN de mai jos /* HIDLIN: Un program simplu pentru a elimina liniile ascunse */ /* Rezultatul acestui program este scris într-un fișier */ /* A ZCRAT pentru a fi citit */ /* Programul GENPLOT */ /* Un program simplu pentru hldden-ІІнѳ ellmlnatlon */ /* Ieșirea acestui program este fișierul A SCRATCH, */ /* care urmează să fie citit de GENPLOT */ tflnclude "stdlo h" #include "math h" tfdefline NVERTEX #define NTRIANGLE unslgned Int STACK- ; /* Latice C */ int ntr-O; dublu v , v , v v v , v v , v , v , eps- e- , meps e- , oneminus- - e~ , oneplus- + e- ; FIȘIER *fout; struct { float X Y; intcode; }s; Capitolul ELIMINAREA LINIILOR ASCUNSE struct { dublu x, y z; } VERTEX[NVERTEX]; struct { int A B, C; dublu a, b, c, h; } TRIUNG[NTRIUNGHI]; maln(argc, argv) Int argc; char *argv[ ]; { Intl, ABCP Q; dublu xO, yO, zO, rho, theta, phl, x, y, z, xA, yA, zA, xB, yB, zB, xC, yC, zC, a, b, c, h, r; FIȘIER *fpin; inch; Dacă (argc!- II (fpln-fopen(argv[ ], "r")) NULL) errmess /* "Fișierul de intrare nu este specificat corect" */ ("Fișierul de intrare nu este specificat corect"); fscanf(fpin, "%lf %lf %lf", &xO, &yO, &zO); skipf(fpln); fscanf(fpin, "%lf %lf %lf", &rho, &theta, &phl); coeff(rho, theta , phi); în timp ce (sklpf(fpln), ch-getc(fpln), ch!-'#' && ch!-EOF) { ungetc(ch, fpln); fscanf(fpin, "%d %lf %lf %lf", &l, &x, &y, &z); If (KO II i>-NVERTEX) errmess /* "număr de vârf ilegal"*/ ("Număr de vârf invalid"); vizualizare(x-xO, y-yO, z-zO, &VERTEX[i] x, &VERTEX[l] y, &VERTEX[l] z); Dacă (VERTEX[i] z ) { if (ntr r NTRIANGLE) errmess /* "Prea multe triunghiuri" */ ("Prea multe triunghiuri"); g - sqrt(a*a+b*b+c*c); a - a / r; b - b / r; c - c/r; h-h/r; TRIANGUL[ntr] A - A; TRIANGUL[ntr] B - B; TRIANGUL[ntr] C - C; ALGORITM PENTRU DETERMINAREA LINIILOR INVISIBILE TRIUNGH[ntr] a - a; TRIANGUL[ntr] bb; TRIUNGH[ntr] c - c; TRIANGUL[ntr++] h - h; } /* Dacă h- , planul ABC trece prin E și hldes nothlng */ /* Dacă h ) llnesegment(VERTEX[P] x, VERTEX[P] y, VERTEX[P] z, VERTEX[ Q] x, VERTEX[Q] y, VERTEX[Q] z, ); fclose(fout); } skipf(fpin) FIȘIER *fpin; {int ch; în timp ce (ch-getc(fpin), ch !- '\n' && ch !- EOF); ț errmess(str) char *str; { printf("%s\n", str); exlt( ); } coeff(rho, theta, phl) dublu rho, theta, phi; { dublu th, ph, costh, slnth, cosph, slnph, factor; factor-atan( , )/ , ; /* Unghiuri în radiani: */ /* Unghiuri în radiani */ factor-th-theta*; ph-phi*factor; cost-cos(th); sinth-sin(th); cosph-cos(ph); sinph-sln(ph); v -sintet; v -cost; /* Elementele matrlx V, vezi Eq ( ): */ /* Elemente ale matricei V, vezi ( ): */ v -cosph*costh; v -slnph*cost; v -cosph*sinth; v -sinph*slnth; v -sinph; v -cosph; v -rho; } vizualizare(x, y, z, pxe, pye, pze) dublu x, y, z, *pxe, *pye, *pze; { /* Coordonatele ochiului, calculate ca în Ec ( ): */ /* Vedeți coordonatele calculate de ( ): */ Capitolul ELIMINAREA LINIILOR ASCUNSE *rhe - ѵ *x + ѵ *y; *rue - v *x + v *y + v *z; *pze - v *x + v *y + v *z + v ; } Hnesegment(xP, yP, zP, xQ, yQ, zQ, JO) xP dublu, yP, zP, xQ, yQ, zQ; intJO; { /* Segmentul de linie PQ Urmează să fie desenat atat cat este */ /* nu este susținut de tranglele JO, jO+ ntr */ /* Segmentul de linie PQ trebuie trasat deoarece */ /* nu este închis de triunghiuri JO JO+ ntr */ Int j-JO, worktodo- , A, B, C, I, Pbeyond, Qbeyond, outslde, Poutslde, Qoutslde, eA eB, eC, suma; dublu a, b, c, h, hP, hQ, r , r , r , xA, yA, zA, xB, yB, zB, xC, yC, zC, dA, dB, dC, MIN MAX, laborator, mu xmin, ymln, zmln, xmax, ymax zmax, C , C , C , K K Scurtcircuit, denoml, denom , Cpos, Ppos, Qpos, aux, epsl; whlle(j eps ? : dA eps ? : dB eps ? : dC - ) {j++; continua; } ALGORITM PENTRU DETERMINAREA LINIILOR INVISIBILE /* Dacă acest test reușește, linia (infinită) PQ lles */ /* în afara piramidei EABC (cel mult un punct comun) */ /* Dacă testul eșuează, există un punct de Intersectlon */ /* Dacă acest test reușește, atunci (infinit) */ /* linia dreaptă PQ se află în afara piramidei EABC */ /* (sau are cel puțin un punct în comun) */ /* Dacă testul eșuează, atunci există un punct de intersecție */ /*Te$t */ /*Test */ Poutslde-Qoutslde-Yu; MIN- ; MAHCH) ; pentru (I-O; KZ; I-n-) { C -yA*zB-yB*zA; C -zA*xB~zB*xA; SZ-xA*uV-xV*uA; /*C x + C y + C z- ls plan EAB */ /* ecuația descrie planul EAB */ Cpos-C *xC+C *yC+C *zC; Ppos-C *xP+C *yP+C *zP; Qpos-C *xQ+C *yQ+C *zQ; denumirea -Qpos-Ppos; If (Cpos>eps) { Pbeyond- Ppos eps; Qbeyond- Qpos>eps; outslde- Pbeyond && Qpos>-meps Qdincolo de && Ppos>-meps; } else outslde- ; Dacă (în afară) se rupe; Iab-fabs(denom ) -mepș && mu -meps &&lab MAX) MAX-lab; } aux-xA; xA-xB; xB-xC; xC-aux; aux-yA; uA-uV; uV-uS; yC-aux; aux-zA; zA-zB; zB-zC; zC-aux; aux-dĂ; dA-dB; dB-dC; dc-aux; } Capitolul ELIMINAREA LINIILOR ASCUNSE Dacă (în afara) {j++; continua;} /*Test */ /*Test */ Dacă (!(Poutslde II Qoutslde)) { worktodo-O; pauză; /* PQ Invlslble */ /* Segmentul PQ invizibil */ } /*Test */ /*Test */ r -xQ-xP; r -yQ~yP; r -zQ-zP; xmin-xP+MIN*r ; ymin-yP+MIN*r ; zmln-zP+MIN*r ; Dacă (a*xmln+b*ymln+c*zmln 'în X y O POLIGONI ȘI PIXELI Dacă punctele C și V dau aceleași semne cu această înlocuire, atunci ele sunt de aceeași parte a liniei A, ceea ce înseamnă că această latură este deasupra Intrarea C este foarte scurtă: topcode[ ] = (D * DAB > ); unde ξ>' și DAB denotă determinanții HA Wah o \u d ■ *: în SW XC Noi HA Wah °AB \u d hv Uv m Determinanții se calculează în mod similar HA Wah dac~ m XC Noi m dbc = xb ^b xs noi după care obținem elementele rămase ale matricei topcode[l]= (D *SHS> ); topcode[ ]= (D *DBC> ); Următorul pas este de a determina coordonatele punctelor finale din stânga (Xleft[l ], Yleft[ ]) și din dreapta (Xright[l ], Yright[l ]) pentru fiecare parte a lui I (/= , , ) Deci, pentru fig vom avea Xleft [ ] = XB Yleft[ ] = YB și așa mai departe Capitolul ÎNCĂRTAREA LINIILOR ASCUNSE Să găsim numerele ipxtype și ipxmax ale coloanelor de pixeli, în care se află cele mai extreme vârfuri din stânga și dreapta ale triunghiului ABC După aceea, pentru triunghiul ABC, puteți ajusta lista de pixeli, limitând-o la numai pixeli din intervalul de numere de coloane de la ipxtype la ipxmax Pentru fiecare coloană ipix din acest interval, introducem margini de-a lungul liniilor LOWER[ipix] și UPPER[ipix] Toți pixelii din intervalul LOWER[ipix], , UPPER[ipix] sunt asociați triunghiului Acest interval include subintervalul (posibil gol) LOW[ipix]+ , , UP[ipix}-i Toți pixelii din această subgamă sunt acoperiți complet de triunghi Pe fig arată latura numărul a triunghiului ABC Să presupunem că această latură este limita inferioară, adică pentru topcode[ \ ] la , atunci această latură a triunghiului vă va permite să definiți valorile LOWER {ipix ] și LOW[ipix ] pentru valorile ipix din interval ipixleft, , ipixright, unde numerele întregi ipixleft și ipixright rezultă din trunchierea coeficientului la calcul ipixleft = (Xleft [l ] - Хтіп) / deltaX \ ipixrigt "(Xright [Z ] - Xmin) IdeltaX', Calculați panta laturilor triunghiului pantă = (Yright [Z ] - Yleft [Z ]) / (Xright [Z ] - Xstânga [Z ]); ' (Spre deosebire de cele două calcule anterioare, aici nu va exista nicio tăiere, deoarece panta variabilă este de tip virgulă mobilă ) Să presupunem că evaluarea secvențială de la coloana ipixleft la coloana ipixright ia în considerare coloana ipx Să determinăm numărul dreptei jl pentru punctul I din fig Orez Latura triunghiului și coloane de pixeli POLIGONI ȘI PIXELI , în care latura triunghiului intersectează limita dintre coloanele de pixeli ipіx și ipіx + Acest număr de linie poate fi găsit din coordonatele ecranului (XI YD ale I: ХІ-Хтіп + (Фіх+ ) * deltaX* YI = Ystânga[l] +pantă* (XI-Xstânga[l]); jl = (YI- Ymin)/deltaY-, /* trunchiat implicit */ /* tăiere implicită */ Fie j old valoarea veche // calculată pentru numărul coloanei ipix - Atunci avem LOWER [ipix] = min(j vechi, jl) ; LOW[ipix ] - max(j old, jl); unde funcția min determină valoarea minimă a unuia dintre cele două argumente, iar funcția max determină valoarea maximă Pe fig avem LOWER[ipix] = , LOW[ipix] = Pentru coloana pixeli ipixleft atribuim coeficientul trunchiat (Yleft[/] - Ymin)/ deltaY iar dacă este numărul liniei pentru punctul final din stânga, atunci jjold În mod similar, numărul de linie trunchiat (Yright[l] - Ymin) /deltaY al punctului final din dreapta este luat în coloana din extrema dreaptă, ipixright Metoda descrisă nu dă întotdeauna rezultatul corect Deci, dacă în fig pentru a considera laturile si in aceasta ordine, atunci obtinem un rezultat incorect: LOWER[ipix] = , deoarece valoarea exacta a lui obtinuta la analiza laturii se va pierde Orez Valoarea LOWER[ipix] depinde de cele două laturi ale triunghiului Capitolul ÎNDEPARTAREA LINIILOR ASCUNSE apoi rezultatul analizei părții Prin urmare, la început vom atribui valori foarte mari tuturor elementelor matricei LOWER și vom permite doar să le reducem De asemenea, tuturor elementelor matricei LOW trebuie să li se atribuie valori foarte mici, care pot fi apoi crescute Valorile matricei UPPER nUP sunt definite în mod similar Acum pentru triunghiul ABC (ale cărui date sunt înregistrate în tabloul TRIANGLE[j]), intervalele coloanelor ipixtype, , ipixmax și elementele tablourilor LOWER[ipix], UPPER [ipix], LOW [ipix], UP [ipix] (ipixtype tfinclude tfinclude tfdefine max (x,y) (W>(y)?(x):(y)) tfdefine min (x,y) ((x) (y) ? max (x,z): max (y,z)) #define min (x,y,z) ((x) tr cov- ; pointer->tr dist-big; po I n te r->s ta rt-N U LL; } reflo(&xO); reflo(&yO); reflo(&zO); printf /* Dați coordonatele sferice rho, theta, phl ale */ /* vlewpolnt E (unghiul phl dintre axele z și OE) */ ("Setați coordonatele sferice rho, theta, phl"); printf("pentru punctul de observatie E\n(phl - unghi intre"); prlntf("axa z si directia segmentului OE):\n"); scanf("%lf %lf %Jf" &rho &theta, &phl); coeff(rho, theta, phl); lnlt vlewport(&Xvp mln, &Xvp max, &Yvp mln, &Yvp max); /* Initlallze vertex array */ /* Init vertex array */ for(l- ; Knvertex; IH-) VERTEX[l] connect-NULL; /* Citește verticile */ /* Citește coordonatele vârfurilor */ xmln-ymin-blg; Xmax-Ymax-mare; whlle (sklpbl(), ch-getc(fpin), ch!-'F' && ch!-'f') { ungetc(ch, fpin); relnt(&l); reflo(&x); reflo(&y); reflo(&z); If (KO II l>-nvertex) errmes /* illegal vertex number */ ("Numărul vârfurilor este negativ sau mai mare decât maximul ( )"); vlewlng(x-xO, y-yO, z-zO, &xe, &ye, &ze); Daca (ze Xmax) Xmax-X; dacă (Y Ymax) Ymax-Y; VERTEX[l] x-xe; VERTEX[i] y-ye; VERTEX[i] z-ze; VERTEX[i] connect - ptr - (int*)malloc(isize); if (ptr* -NULL) errmes /* "Eroare de alocare memorie " */ ("Eroare de alocare memorie "); *ptr ; } /* Calculează constantele ecranului */ /* Calculează constantele ecranului */ Xrange-Xmax-Xmin; Yrange-Ymax-Ymin; Gama Xvp-Xvp max-Xvp min; Gama Yvp-Yvp max-Yvp min; fx-Xvp range/Xrange; fy-Yvp range/Yrange; PROGRAM ÎMBUNĂTĂTIT d-(fx ) {POLY[ }-i; npoli- ; sklpbl(); în timp ce (ch-getc(fpin), ch !- '#') { ungetc(ch, fpin); reint(&POLY[npoli-H-]); if (npoly NPOLY) errmes /* "Prea multe vârfuri într-un singur poligon" */ ("Prea multe vârfuri într-un singur poligon"); } If (npoly- - ) errmes /* "Only one vertex of polygon" */ ("Un poligon are un singur vârf"); Dacă(npoli ) { add linesegment(POLY[ ], POLY[ ]); continua; } dacă (!counter clock( , , , &dlag, )) continuă; /* Fața din spate */ /* Fața din spate */ pentru ( - ; i ) { min dlag-big; pentru( -Yu; IKnpoly; ++) { IO-(I ? npoli- : - ); I -( npoli- ? : + ); ** Capitolul ELIMINAREA LINIILOR ASCUNSE if (contor ceas(IO AND, i , &diag,O) && diag - ) { pnode-^struct node *)malloc(sizeof(struct node)); if (pnode NULL) errmes /* "Eroare de alocare memorie " */ ("Eroare de alocare memorie "); pnode->jtr-pointer~>tr cov; pnode->next-polnter->start; pointer->start-pnode; } } /* Desenați toate segmentele de linie cât de mult sunt vizibile */ /* Desenați segmente de linie dreaptă când sunt vizibile */ pentru (P- ; P conectare; if (ptr - - NULL) continua; xP-pvertex->x; yP - pvertex->y; zP - pvertex->z; XP-xP/zP; YP-yP/zP; pentru (lconnect- ; lconnect x; yQ - pvertex->y; zQ - pvertex->z; XQ-xQ/zQ; YQ-yQ/zQ; /* Folosind listele de ecran, vom extinde */ /* set de triunghiuri care pot conține polenți de PQ: */ /* Pe baza listelor de ecran, definiți un set de triunghiuri */ /* care poate acoperi punctele segmentului de dreaptă PQ: */ PROGRAM ÎMBUNĂTĂTIT Dacă (XP start; whlle (pnod!-NULL) { trnr* pnode->jtr; /* trnr va fi stocat doar dacă nu este încă */ /* prezent în matrice trset (setul trangle) */ /* Numărul trnr va fi reținut numai dacă */ /* nu este deja în tabloul trset (un set de triunghiuri) */ trset[ntrsetpirnr; /* santinel */ /* "santinela" */ jtrO; whlle (trset[|tr]!-trnr)jtr++; Dacă(jtr ntrset) { ntrset-n-; /* asta înseamnă că trnr este stocat */ /* asta înseamnă că trnr este scris într-o matrice */ Dacă (ntrset- -nntrset) errmes /* "Trangle set overflow" */ ("Setul de triunghiuri este plin\n"); } pnode-pnode->next; } } /* Acum trsetf ] trset[ntrset- ] Este setul de */ /* trangles care pot hlde polnts de PQ */ /* Acum lista trsetfO] trset[ntrset- ] definește */ /* set de triunghiuri care pot acoperi punctele PQ */ llnesegment(xP, yP, zP, xQ, yQ, zQ, ); } } endgrf); } Capitolul ELIMINAREA LINIILOR ASCUNSE sklpbl() { char ch; do ch-getc(fpin); în timp ce (Isspace(ch) l comment(ch)); ungetc(ch fpln); } /* */ int comment(ch) char ch; { intk; Dacă (ch 'C) { do k-getc(fpln); în timp ce (k !- ')' && la !- EOF); returnează k- -')'; } altfel returnează ; } /* */ int reflo(px) dublu *px; { sklpbl(); returnează fscanf(fpln, "%lf", px); } /* */ int reint(pi) int *pi; { skipbl(); returnează fscanf(fpin, "%d", pi); } /* */ -add Hnesegment(P, Q) Int P, Q; { int iaux, *ptr, ii, n; Dacă (P>Q){ iaux-P; PQ; Q-lux } /*Acum:P : verificați dacă următorul triunghi este coplanar; depoziteaza-l */ /* */ /* cod - ; orientare: dacă în sens invers acelor de ceasornic - */ /★ calculul lungimii proiecției diagonalei AC */ PROGRAM ÎMBUNĂTĂTIT /* cod - : calculează a b, c, h; memorarea primei */ /* triunghiuri */ /* cod > : verificați coplanaritatea următoarelor */ /* triunghi; scrie-l în memorie */ { Int A-abs(POLY[IO]) B-abs(POLY[i ]) C-abs(POLY[l ]); dublu xA, yA, zA, xB, yB, zB, xC, yC, zC, r, xdist, ydist, zdlst XA Y A XB YB XC Y C hO D A D B DC D DAB, DAC DBC, aux, dlst, xR yR; static dublu abc h; pvertex-VERTEX+A; xA-pvertex->x; yA - pvertex->y; zA - pvertex->z; pvertex-VERTEX+B; xB - pvertex->x; yB - pvertex->y; zB - pvertex->z; pvertex-VERTEX+C; xC-pvertex->x; yC - pvertex->y; zC - pvertex->z; hO-xA* (yB*zC-yC*zB)-xB* (yA*zC - yC*zA)+ xC* (yA*zB - yB*zA); Dacă (cod- -Yu) Dacă (ho>eps) { xdlst-xC-xA; ydlst-yC-yA; zdlst-zC-zA; *pdist-xdlst*xdist+ydist*ydist+zdist*zdist; întoarcere ; } altfel returnează ; /* Dacă h - , planul ABC trece prin E și nu conține nimic */ /* Dacă h A - A; triunghi->B - B ; triunghi->C - C; Capitolul ELIMINAREA LINIILOR ASCUNSE triunghi->a - a; triunghi->b - b ; triunghi->c - c; triunghi->h - h; /* Triunghiul va fi acum stocat în ecranul llsts */ /* din pixelii asociați; mai întâi tablourile LOWER, */ /* SUS SCĂZUT UP sunt definite: */ /* Acum triunghiul trebuie să fie scris în listele de ecran */ /* de pixeli asociați, dar trebuie făcut mai întâi */ /* definește matrice LOWER, UPPER, LOW, UP: */ XA-xA/zA; YA-yA/zA; XB-xB/zB; YB-yB/zB; XC-xC/zC; YC-yC/zC; DA-XB*YC-XC*YB; DB-XC*YA-XA*YC; DC-XA*YB-XB*YA; D-DA+DB+DC; DAB-DC-M*(XA-XB); DAC-DB-M*(XC-XA); DBC-DA-M*(XB-XC); topcode[ ]-(D*DAB> ); topcode[ ]-(D*DAC> ); topcode[ ]-(D*DBC> ); Xstânga[ }-XA; Yleft[ }-YA; Xdreapta[O}-XB; Yright[O}-YB; Xleft[ ]-XA; Yleft[ }-YA; Xdreapta[ }-XC; Yright[ }-YC; Xstânga[ }-XB; Ystânga[ }-YB; Xrlght[ }-XC; Yright[ }-YC; formă); l Xright[l] (XleftflJ- -Xrightfl] && Yleft[l>Yright[l])) { aux-Xleft[l]; Xstânga[l}-Xdreapta[l]; Xright[l}-aux; aux-Yleft[l]; Ylefl(l}-Yright[l]; Ylefl[l}-aux; } ipixmin-xwhole(min (XA,XB,XC)); ipIxmax-xwholeJmaxâfXA XB XC)); pentru (ipix-iplxmln; ipix LOW[ipix] && jplx eps ? h*sqrt(xR*xR+yR*yR+ )/denom : mare; /* Linia de la punctul vlew E la plxel polnt(xR, yR, ) */ /* intersectează planul ABC la o distanță dlst față de E */ /* Linia de la punctul ochiului E la punctul pixelului (xR, yR, ) */ /* intersectează planul ABC la distanța dlst de E */ Dacă (dlst tr dlst) { pointer->tr cov-ntr, pointer->tr dist-dlst; } } else /* Adăugați triunghiul pe ecran llst: */ /* Triunghiul este adăugat la lista ecranului */ { pnode-^struct node *) malloc(sizeof(struct node)); If (pnode- -NULL) errmes /* "Eroare de alocare memorie " */ ("Eroare de alocare memorie "); pnode->jtr - ntr; pnode->next - polnter->start; pointer->start - pnode ; } } ntr-n-; errmes(str) char *str; {fatal(); printf("%s\n", str); exlt( ); /* fatal Este deflned În același fișier ca endgr etc */ /* fatal este definit în același fișier ca endgr etc */ } /* -*/ coeff(rho, theta, phl) dublu rho, theta, phl; { dublu th, ph, costh, slnth, cosph, slnph, factor; factor-atan( , )/ , ; /* Unghiuri în radlani */ /* Unghiuri în radiani */ th-theta*factor; ph-phl*factor; cost-cos(th); slnth-sin(th); Capitolul ELIMINAREA LINIILOR ASCUNSE cosph-cos(ph); sinph-sln(ph); /* Elementele matricei V, vezi Eq ( , ): */ /* Elementele matricei V, vezi ( ): */ v -slnth; v -cosph*cost; v -sinph*cost; v -cost; v -cosph*sinth; v ^>-sinph*slnth; v -sinph; v -cosph; v -rho; vlew ng(x, y, z, pxe, pye, pze) dublu x, y, z, *pxe, *pye, *pze; { /* Coordonatele ochiului, calculate ca în Ec ( ): */ /* Vedeți coordonatele calculate de ( ): */ *pxe-v *x +v *y; *rue - v *x + v *y + v *z; *pze - v *x + v *y + v *z + v ; } /* */ llnesegment(xP, yP, zP, xQ, yQ, zQ kO) dublu xP, yP, zP, xQ, yQ, zQ; Int kO; { /* Segmentul de linie PQ trebuie trasat, în măsura în care nu este */ /* susținut de triunghiuri trset[kO] la trset[ntrset- ] */ /* Segmentul de linie dreaptă PQ este trasat deoarece */ /* nu este închis de triunghiuri de la trsetfkO] la trset[ntrset- ] */ Int j, k-kO, worktodo- , A, B, C, I, Pbeyond, Qbeyond, outslde, Poutside, Qoutside, eA, eB, eC, sum; dublu a, b, c, h, hP, hQ, rl, r , r , xA, yA, zA, xB, yB, zB, xC, yC zC, dA, dB, dC, labmln, labmax, lab, mu, xmin, ymin, zmin, xmax, ymax, zmax, C C C K K KZ denoml, denom , Cpos, Ppos, Qpos, aux, eps ; în timp ce (k a; b-triunghi->b; c-triunghi~>c; h-triunghi->h; /* Testul */ /* Testul */ hP-a*xP+b*yP+c*zP; hQ-a*xQ+b*yQ+c*zQ; eps -eps+eps*h; dacă (hP-h A; B-triunghi->B; C-ptrangle->C; pvertex-VE RTE X+A; PROGRAM ÎMBUNĂTĂTIT xA-pvertex->x; yA - pvertex->y; zA - pvertex->z; pvertex-VERTEX+B; xB - pvertex->x; yB - pvertex->y; zB - pvertex->z; pvertex-VERTEX+C; xC-pvertex->x; yC - pvertex->y; zC - pvertex->z; dA-K *xA+K *yA+K *zA; dB-K *xB+K *yB+K *zB; dC-K *xC+K *yC+K *zC; /* Dacă dA, dB, dC au același semn, vârfurile */ /* A, B, C se află pe aceeași parte a planului EPQ */ /★ Dacă dA, dB, dC au aceleași semne, atunci vârfurile ★/ /★ A, B, C sunt de aceeași parte a planului EPQ */ eA-dA>eps ? : dA eps? : dB - ) { k-n-; continua; } /* Dacă acest test reușește, linia (infinită) PQ se află */ /* în afara pyrfcmid EABC (sau linia și piramida */ /* au cel mult un punct în comun) */ /* Dacă testul eșuează, există un punct de intersectlon */ /* Dacă acest test reușește, atunci (infinit) */ /* linia dreaptă PQ se află în afara piramidei EABC */ /* (sau are cel puțin un punct în comun) */ /* Dacă testul eșuează, atunci există un punct de intersecție */ /*Test */ /*Test */ Poutside-Qoutslde-Yu; labmin- ; Iabmax-Yu ; pentru (i- ; i eps) { Pbeyond- Ppos eps; Qbeyond- Qpos>eps; afară- Pbeyond && Qpos>-meps Qbeyond && Ppos>-meps; } else outside- ; Dacă (în afară) se rupe; lab-fabs(denom ) -meps && mu -meps && lab - ? p+bază : p-bază; } Pentru a înțelege acest program, puteți consulta și paragraful , care descrie fișierul de intrare pentru desenarea unei singure litere L Acum, în loc de aceasta, linii de formă numărul vertex x y z în ordinea crescătoare a numerelor de vârfuri, deși ordinea în care sunt enumerate nu are efect asupra calculului Pentru a evita apelurile multiple la funcția fprintf, au fost introduse funcțiile wl și w Ele, la rândul lor, se referă la funcția s, care oferă o creștere a numerelor inițiale de vârfuri ( -n+ ; i ) fprintf(fp, "%d\n", i); Capitolul EXEMPLE PRACTICE fprintf(fp, "%d %d#\n", *n,-n); /* Față de delimitare inferioară */ /* Față de delimitare inferioară */ fprlntf(fp, "%d %d\n", *n,- *n); for(l- *n+ ; K- *n; I++) fprlntf(fp, "%d\n", I); fprlntf(fp, "%d\n", *n); for(l- *n- ; l>- *n+ ; I ) fprlntf(fp, "%d\n", I); fprlntf(fp, "%d#\n", *n+ ); /* Linii verticale */ /* Linii verticale */ pentru ( - ; Kn; I++) { JI%n+ ; fprlntf(fp, "%d %d %d %d#\n", J, I l+ *n, J+ *n); fprintf(fp, "%d %d %d %d#\n", l+n, j+n, j+ *n, l+ *n); } fclose(fp); } TIJE ÎN SPIRALĂ Pe fig arată următorul exemplu - o spirală formată din tije orizontale de lungime Z, lățime și înălțime w Tija inferioară se află pe planul xy, așa cum se arată în Fig Începând de jos, poziția fiecărei tije următoare se obține prin rotirea celei anterioare în jurul axei z cu un unghi de ° și un incrementul de timp al coordonatei sale z cu xv Numerotăm tijele , , , n - de jos în sus Tija i are numerele de vârfuri /, / + , , z + atribuite ciclic, în conformitate cu numerele de vârfuri ale tijei prezentate în Fig Fiecare punct (l', y', z') al barei i + poate fi obţinut prin rotirea punctului corespunzător (x, y, z) al barei i cu ° în jurul axei z (şi stabilirea z* = z + w), atunci există [x' y']= [x y] cos ° -sin ° sin ° cos ° sau doar x' = -y și y' = x Acesta este ceea ce este folosit în următorul program /* Preprocesor BEAMS pentru HIDLINPX) */ /* preprocesor pentru programul HIDLINPX */ #include "stdio h" maln() { FIȘIER *fp; Int I J, n, A; float I, w, xA, yA, xB, yB, xC, yC, xD, yD, a, b, aux, z; printf /* "Câte grinzi?" */ ("Câte bare?\n"); scanf("%d", &n); Capitolul EXEMPLE PRACTICE printf /* "Fascicul măsoară I xwx wAnGive I și w:" */ ("Fascicul măsoară Ixwx ѵѵAnDați valorile lui I și w:"); scanf("%f %f", &l &w); fp-fopen("beams dat", "w"); fprlntf(fp, " , , %f\n", n*w/ , ); /* obiect central polnt */ /* punct central al obiectului ★/ a- , *I; baw; xA-a; uA~-a; xV-a; UV-a; xC-b; yC-a; xD-b; yD-a; pentru (i- ; i r) Pașii leagă stâlpul de tijele verticale Etapa de jos este prezentată în fig Fiecare treaptă este un bloc de lungime R~r, lățime , A și înălțime , L Împreună, toți pașii fac o revoluție completă, astfel încât unghiul de rotație pentru fiecare pas este l P Capitolul EXEMPLE PRACTICE Vom atribui pașilor numerele , , n - , numerotându-le de jos în sus Pasul i poate fi obținut prin rotirea treptei zero în jurul axei z cu un unghi a = id și ridicând-o simultan la o înălțime ih Astfel, coordonatele vârfurilor (x, y, z) ale pasului i pot fi calculate din coordonatele vârfurilor corespunzătoare (X, Y, Z) pasului : [xy]-[X Y] cosa șina -șina cosa z = Z + ih În programul de mai jos, coordonatele vârfurilor treptei inferioare sunt scrise în rețele X, Y, Z Vârfurilor treptelor, inclusiv tijele verticale și balustradele, li se atribuie numerele , , p - Presupunând L / \u d p, puncte echidistante ale cercului inferior al coloanei centrale, atribuim următoarele n numere lui M, , M + n - În mod similar, numerele întregi M + n, , M + - sunt atribuite cercului superior al coloanei SCARA IN SPIRAOLA /* SĂRĂ-ȘIRURĂ (preprocesor pentru HIDLINPX) */ /* Preprocesor pentru programul HIDLIPX*/ ♦Includeți "stdlo h" ♦include "math h" principal() { FIȘIER *fp; Int I, J, n, k, M; float r, R, pi, alpha, cosa, șina, x, y, z, delta, h, H, X[ ], V[ ] Z[ ]; printf /* "Dă numărul n (pentru a desena n stalri)" */ ("Specificați numărul n (pentru a desena n pași):\n"); scanf("%d", &n); printf /* "Dă ajutor unei singure trepte de scară" */ ("Setați înălțimea unei trepte individuale:"); scanf("%f", &h); printf /* "Dați raza mare R și raza mică r\n" */ ("Precizați razele mari (R) și mici (r) ale scărilor\n"); scanf("%f %f", R, &r); fp-fopen("wlndlng dat", "w"); pl- *atan( ); delta- *pi/n; fprlntf(fp, " , , %f\n", , *n*h); pentru(j- ;j h/ ; X[ ]-R; Y[ > ; Z[ }- *h; M- *n; Hn*h; pentru H); l -M; i ) fprintf(fp, "%d", I); fprlntf(fp, "#\n"); /* partea de jos a câmpului */ /* partea de jos a coloanei */ for(i-M+n; l r) Pe fig un cerc orizontal mare determină poziția centrelor cercurilor care formează torul, raza acestui cerc este R Să alegem n puncte echidistante de pe acest cerc ca centre de verticală mică Orez Thor TOP cercuri situate cu raza r Reprezentarea parametrică a unui cerc mare este descrisă prin formule x = R cos a y = Rșina z = Punctul corespunzător lui a = este centrul cercului mic x \u d R + r cos fi y \u d z \u d z sin P care este prezentată și în Fig Restul n - cercuri mici se formează prin rotirea acestui cerc inițial în jurul axei z cu unghi = d, unde i= , ,m- u și y= sunt numerotate , , , m, numărând de jos în sus Punctele dintr-un semicerc vecin pentru y > sunt numerotate n + , n + , , m și așa mai departe Acest mod de numerotare este ilustrat în tabelul următor, unde fiecare linie (orizontală) corespunde punctelor de pe orizontală Orez emisferă Capitolul EXEMPLE PRACTICE cerc și, în mod similar, fiecare coloană corespunde punctelor de pe un sfert din cerc vertical: i -* m- / p p p Zm m m + m + ( m- )m + m + m+ ( m- )m + Stratul inferior al emisferei este format din triunghiuri, toate având un vârf comun în punctul numărul Restul n - straturi constau din patrulatere ABCD, fiecare dintre ele având două margini orizontale paralele reciproc AB și DC Rețineți că această emisferă este fundamental diferită de toate celelalte exemple considerate mai devreme Acesta nu este un obiect solid, ci o suprafață vizibilă din orice direcție, în funcție de poziția punctului de observare Deoarece fiecare dintre patrulaturi are două laturi, programul trebuie să descrie ambele secvențe, ABCD și DCBA /* SEMI SPHERE (preprocesor pentru HIDLINPX) */ /* preprocesor pentru programul HIDLINPX */ #include "stdio h" #include "math h" main() { FILE *fp; Int i, j n, A, B, C, D, P Q; float pi, alpha, beta, delta, cosa, șina, cosb, sinb; printf /* "Dă numărul n" */ ("Specificați numărul n\n"); scanf("%d", &n); fp-fopen("semi dat", "w"); pi- , *atan( , ); delta-pi/( *n); /* n * delta - pl/ */ fprintf(fp, " - \n"); /* punct central obiect */ /* punctul central al obiectului */ /* R - ; centrul sferei în O */ /* Raza - unitate, centrul sferei la O */ fprintf(fp,"O - \n"); /* primul punct */ /* primul punct */ FUNCȚIA A DOUĂ VARIABILE pentru H); K *n; I++) {alfa-l*delta; cosa-cos(alfa); sina-sln(alfa); pentru - ; j = , ) {/= ; / = ;} altfel {R = ; /= ;} dacă {u = , ; apoi operatorul iv = , ; ar fi întotdeauna executată, deoarece această instrucțiune nu are nimic de-a face cu condiția și MAX) mergeți la gata, z++;} gata același efect poate fi obținut fără a utiliza declarația necondiționată goto, de exemplu i = ; în timp ce (i MAX) pauză*, z++;} sau pentru (z = ; i MAX) pauză ;} Puteți utiliza operatorul de continuare pentru a trece imediat la verificarea condiției de terminare a buclei continua; De exemplu, operatorul în timp ce (a Alocarea unui element al structurii adresat de pointer (Următoarele operații sunt de obicei numite unare: au un singur operand Rețineți că simbolurile ♦ sunt folosite și ca operații binare radiouri) Negație logică DREAPTA LA STÂNGA negație pe biți ++ Schimbarea semnului Crește cu unu - Scade cu unu - & Definirea adresei * Contact prin adresa (tip) Type conversion sizeof Determinați dimensiunea în octeți * Înmulțire de la stânga la dreapta / Divizia % Modulo (restul) + Adăugare de la stânga la dreapta - Scăderea " Shift Stânga Stânga Dreapta " Schimbați la dreapta mai mult decât >= Mai mare sau egal Aplicație O SCURTĂ INTRODUCERE LA C != Egal Nu este egal De la stânga la dreapta & Operare pe biți AND De la stânga la dreapta Operațiune XOR De la stânga la dreapta Operare pe biți SAU De la stânga la dreapta && Operațiune logică AND De la stânga la dreapta II Operare logică SAU De la stânga la dreapta Operațiune condiționată DREAPTA LA STÂNGA Atribuire DREAPTA LA STÂNGA , Operațiune virgulă De la stânga la dreapta O expresie nu numai că poate schimba valori, ci și poate efectua acțiuni care schimbă starea Următoarea linie, de exemplu, este expresia i = i + Aceasta este o expresie de atribuire care crește valoarea lui / cu Mai puțin evident, această expresie are și o valoare, și anume noua valoare a lui i Prin urmare, are sens să scrii І= * (z=z+l) + ; care, în special, poate fi înlocuit cu j= * (++i) + ; Ambele expresii ++z și /++ înseamnă că valoarea variabilei este mărită cu unu Cu toate acestea, expresiile produc semnificații diferite: ++z înseamnă variabila incrementală i și utilizare definirea noului său sens; z++ înseamnă folosiți valoarea veche a re- variabila i și apoi incrementând acea variabilă Deci, după efectuarea operațiunilor z = ; /=(++/); t = ; n=(m++); vom avea i=j = , m = , n = A OPERATORI ŞI EXPRIMI O expresie de atribuire și un punct și virgulă (;), în această ordine, formează un operator de atribuire Asa de, / = /+ ; denotă nu o expresie, ci un operator (notația /++; capătă un sens similar) Dacă doriți să efectuați mai multe acțiuni într-un context care permite doar expresii (nu operatori), puteți utiliza operatorul virgulă Luați în considerare, de exemplu, ciclul în timp ce (i + = j, j -, k = i + * j, k> ) j *= ; Funcționează similar cu linia din nou : i + = j ; j-; k = i + * j ; dacă (k> ) {j *= ; du-te din nou ;} (scrierea j *= înseamnă j=j * , j- = înseamnă j=j- și așa mai departe) Câteva personaje? și : formează un alt operator, nu chiar obișnuit, dar foarte convenabil În așa-numita expresie condiționată cond exprl : expr , se evaluează prima expresie cond, ceea ce înseamnă condiția Atunci fie exprl, fie expr sunt evaluate în funcție de faptul dacă valoarea lui cond este diferită de zero sau, respectiv, zero Reamintim că în C, valoarea lui adevărat ("adevărat") este notată cu , iar fals ("fals") cu , deci în loc de dacă (x > Y) max = x; altfel max = y; poate fi scris max = (x > y ? x: y); Este foarte important să se facă distincția între operațiile logice și pe biți Operațiile logice sunt folosite foarte des, de exemplu în această expresie (x && y ) și (y && j/i > k, poate apărea împărțirea la zero Anexă O SCURTĂ INTRODUCERE LA C Rezultatul evaluării expresiilor booleene va fi întotdeauna sau Dar în operațiile pe biți, întreaga secvență de biți a operanzilor este întotdeauna înțeleasă, chiar dacă sunt specificate variabile întregi, cum ar fi i = ; ji" ; k = i\j\ Aici, sunt utilizate operațiile "deplasare la stânga" (") și "operație SAU pe biți" (I) Valorile variabilelor i, /, k pot fi reprezentate în format binar /" (= ) y = (= ) * = (= ) O operație de conversie de tip poate fi utilizată pentru a forța o conversie de tip a unei variabile Deci, sensul expresiei (int) , va fi egal cu numere întregi ("întregi") Dacă ambele variabile і și j sunt de tip întreg, atunci rezultatul împărțirii і /j va fi valoarea întreagă a coeficientului trunchiat Astfel, / are o valoare întreagă , chiar și într-un context în care se așteaptă un float, ca în expresia x = / , unde x este de tip ceremonial real Dacă cel puțin unul dintre operanzii a și b este flotant, atunci câtul a / b este, de asemenea, flotant și nu este trunchiat Constantele în virgulă mobilă conțin o virgulă zecimală sau litera e (sau E); constantele întregi nu le conțin Asa de, e / are valoarea / are o valoare de , / are valoarea Operatorul modulo (%) poate fi aplicat numai operanzilor întregi: % contează % invalid (int) , % are valoarea În ultimul exemplu, operatorul de conversie de tip (int) are prioritate față de operatorul modulo % Dacă o valoare în virgulă mobilă este atribuită unei variabile întregi, va avea loc trunchierea Operator / = , ; A ÎNTREBĂRI LEXICE ȘI STRUCTURA PROGRAMULUI atribuie valoarea variabilei i Dacă variabila x este de tip float, atunci când se execută următoarea instrucțiune, lui x i se va atribui valoarea (transformată într-un număr în virgulă mobilă ) x = / ; A ÎNTREBĂRI LEXICE ȘI STRUCTURA PROGRAMULUI În cele mai multe cazuri, spațiile și liniile noi nu au niciun efect asupra semnificației programului Secvența de caractere din formular /♦ ♦/ este ignorat de compilator Așa sunt indicate comentariile destinate unei persoane Identificatorii, cum ar fi numele de variabile, constau din litere și cifre latine, dar primul caracter trebuie să fie o literă În C, liniuța de subliniere ( ) este de asemenea considerată o literă Ținând cont de acest lucru și făcând distincție între litere mari și mici, numărăm de litere diferite Identificatorii pot fi folosiți și pentru a desemna constante, de exemplu # amendă MAXIM După o astfel de "linie de control al preprocesorului" va fi posibil să se utilizeze identificatorul MAXIM, pur și simplu ca o altă desemnare pentru numărul Acesta este doar primul exemplu de macrocomandă Iată o macrocomandă mai interesantă numită MAX și reprezentată de șir #defineMAX(x, y) x>y? X y ceea ce înseamnă că oriunde în program orice linie a formularului MAX (a, b) înlocuit automat cu șirul a > b? a:b Pentru a evita confuzia în cazuri mai complexe, este recomandat să folosiți paranteze, care pentru această macrocomandă ar putea arăta tdefiniți MAX(x, y) ((x) > (y) ? (x): (y)) Anexă O SCURTĂ INTRODUCERE LA C Există linii de control pentru includerea fișierelor Linie folosită des #include Rezultatul acțiunii sale este înlocuirea acestei linii cu conținutul unui fișier numit stdio h care se numește "fișier antet standard de intrare/ieșire" Dacă sunt utilizate funcții matematice precum cos și sin, atunci o linie de control trebuie inclusă în program #include Programul conține de obicei una sau mai multe funcții Nu există subrutine sau proceduri în limbajul C, ci doar funcții Chiar și programul principal este o funcție numită tip Aceasta este doar o problemă terminologică, deoarece funcțiile nu trebuie să returneze deloc o valoare și pot fi accesate în același mod ca în alte limbi - la o procedură Acest lucru este prezentat în programul următor, care citește coordonatele ortogonale D ale două puncte P și Q și calculează distanța dintre cele două puncte /* Acest program calculează distanța dintre */ /* două puncte date P și Q */ /* Acest program calculează distanța dintre ♦/ /* două date P și Q */ #include principal() {float xP yP xQ yQ; printf ("Setați xP yP xQ yQ:") scanf ('' %/ %/ %/ %G, &xP, &yP, &xQ, &yQ); prinț d ist ance (xP yP xQ, yQ); } print distance (xl, yl, x , y ) float xl, yl, x , y { float delta x, delta y distanţă; delta x = x - xl; delta y = y - yl; distanta = sqrt(delta x * delta x + delta y * delta y); prinț("Distanța: %f\n" distanță); } Funcții standard scanf prinț, scrt va fi descris în paragraful A L MATRII ȘI POINTERS A MATRII ȘI POINTERS După declararea unui tablou //oa/a[ ]; Următoarele cinci variabile float sunt disponibile: a[ ], a[ ], a[ ], a[ ], a[ ] De asemenea, puteți scrie a[r], unde indicele i este orice expresie întreagă a cărei valoare nu poate fi negativă sau mai mare de Indicele încep întotdeauna să conteze de la Într-o declarație de matrice pot apărea numai constante întregi Declarațiile folosesc adesea nume constante în loc de numere, ca în exemplul următor, care arată, de asemenea, că tablourile pot avea mai mult de un index: #define NROWS Zdefiniți NCOLUMNS int tabel[NROWS ] [NCOLUMNS ]; pentru (i = ; i Лп", X [ ] + Y[ ]) } Variabilele statice trebuie să fie precedate de cuvântul cheie static atunci când sunt declarate În următorul exemplu de program, tablourile X și Y sunt interne funcției de tip, dar sunt declarate statice, deci au spațiu constant în memorie și, prin urmare, pot fi inițializate tip() { float static X[ ] = { , , , L}, Y[ ] = {- , , , }; rll/("% Lpi,X[ ]+ Y[ ]); } A FUNCȚII În limbajul C, funcțiile pot face uneori diferența Dacă o funcție trebuie să aibă o valoare returnată, atunci i se atribuie utilizând instrucțiunea return După cum sugerează și numele, apelarea acestui operator face ca funcția să iasă și să iasă imediat Deci următoarele două funcții dau același rezultat int f(x) float x; {dacă(x p unde două caractere - și > seamănă cu o săgeată și formează un singur simbol de operare A ALOCAREA DINAMICĂ A MEMORIEI Să presupunem că a fost declarat char *p, *malloc(), *realloc(); Atunci p va fi un pointer către un simbol și putem scrie p = malloc(n)*, unde n este o expresie întreagă pozitivă care indică numărul de octeți Efectul acestui operator este că, dacă este posibil, o bucată adiacentă de memorie este alocată pentru n caractere Dacă volumul necesar Anexă O SCURTĂ INTRODUCERE LA C memoria nu este disponibilă, atunci variabilei p i se va atribui valoarea NULL - o valoare specială pentru un pointer "null" care nu indică către un obiect real; aceasta vă permite să efectuați o verificare (de exemplu, să emitați un mesaj despre memorie insuficientă) if (p = NULL) { !*Tetogu insuficient */ /*Memorie epuizată */} Cei n octeți alocați în acest fel sunt acum disponibili prin p Să presupunem că acum este necesar să plasăm litera Q în poziția z-a ( n) Pentru a face acest lucru, puteți scrie p = realloctp, N)\ Acum, dacă pointerul p nu este NULL, atunci indică un bloc de N octeți, iar primii n octeți ai acestui bloc vor avea același conținut care le-a fost scris înainte Memoria poate fi alocată și pentru alte tipuri de date Să presupunem că vrem o succesiune de k întregi Deoarece funcția standard malloc trebuie să știe în mod necesar de câți octeți sunt necesari, suntem interesați să știm câți octeți sunt alocați pentru un număr întreg Acest număr este determinat în mod independent de mașină de operație dimensiunea(int) O altă complicație are de-a face cu tipul de indicatori care sunt necesari Avem nevoie de un indicator către un număr întreg în loc de un indicator către un caracter Cu toate acestea, funcția taIos este asociată cu un indicator de caractere Conversia tipului se poate face cu un operator de conversie a tipului Prin urmare, scriem int *p; p = (int *)malloc(k * sizeof(int)); A INTRARE / RETRAGERE Aici operatorul de conversie de tip (zn/*) spune că typedm dorit este un pointer către un număr întreg Numărul întreg j (j= , , ,k - ) din această secvență este notat cu *(p+/) A INTRARE/IEȘIRE Nu există constructe speciale pentru intrare și ieșire în limbajul C În schimb, o serie de funcții standard sunt utilizate împreună cu un tip de structură predefinit numit FILE Toate detaliile despre aceasta sunt incluse în programul nostru într-un mod foarte convenabil folosind linia de preprocesor #include Fișierul antet stdio h conține o declarație typedefb de forma typedef Struct { } FILE\ deci dacă scriem • FILE*fp-, atunci variabila fp va fi de tip "pointer to file FILE", iar computerul va cunoaște toate detaliile despre corpul structurii FILE Structura este de fapt accesibilă prin indicatorul fp după ce fișierul este deschis de către operator fp^fopen(cale-fișier, mod)-, în care argumentele au următoarea semnificație: nume-fișier - un șir care conține numele fișierului așa cum este scris în director; mod ~ caractere "r" sau "w" pentru intrare sau ieșire formatate; pentru I/O brut, formatul de notație depinde de compilatorul utilizat Funcția fopen este simetrică cu funcția fclose-, funcția fopen conectează fișierul la program, iar funcția fclose îl deconectează Următorul exemplu de program arată cum puteți scrie ceva în fișierul EXEMPLU Dacă un astfel de fișier nu a existat înainte, acesta va fi creat #include maini) { FILE'fp-, intz; Aplicație O SCURTĂ INTRODUCERE LA C fp = fopen("EXEMPLU", " ; pentru (/= ; z<= ; /++) fprintf(fp, "i=%ld i*i = % d\n", /, i*i); fclose(fpY, } După executarea acestui program, va fi creat un fișier EXEMPLU cu următorul conținut: Z= z*z= z = z*z = z= z*z= z= z*z= În loc de funcția fprintf, putem defini funcția printf omițând primul argument Apoi rezultatul va apărea pe ecranul de afișare al stației noastre de lucru în loc să fie scris pe disc În acest caz, programul complet ar putea arăta ca tip() { int g, pentru (i = ; z<= ; /++) printfCi =% di*i=Q/Q d\n , i, i*i); } De fapt, printf( ) este echivalent cu fprintf (stdout, ), unde stdout este indicatorul fișierului către ieșirea standard declarată în fișierul antet stdio h Primul argument pentru printf (al doilea argument pentru fprintf) va fi un șir de format care conține porțiunea de text care ar trebui să apară în formă literală în rezultat și elementele de format referitoare la restul argumentelor Aici, elementele de format sunt % d și % Ei definesc că valorile variabilelor i și i*i vor fi numere întregi, care ar trebui tipărite cu una, respectiv două zecimale Toate celelalte caractere din șirul de format sunt tipărite literal, inclusiv spațiile, iar caracterul newline este notat cu \l În elementul de format, litera d trebuie înlocuită cu litera / dacă elementul de date asociat este de tip float în loc de întreg Datele sunt citite atunci când accesați funcții similare fscanf(fp, format-string, ) dacă datele sunt într-un fișier de pe disc, sau scanf(format-string, ) A INTRARE / RETRAGERE dacă datele urmează să fie introduse de la tastatura utilizatorului Cu toate acestea, argumentele rămase trebuie să fie acum pointeri, deoarece funcțiile fscanf și scanf trebuie să poată atribui valori variabilelor prin aceste argumente De exemplu, dacă un număr trebuie specificat de către utilizator, atunci puteți scrie printfС Introduceți un număr: "); scanf C%d\ &număr); (absența simbolului & în desemnarea variabilei ar însemna o greșeală gravă) Dacă variabila x este de tip float și variabila xx este de tip double, atunci valorile acestora se citesc după cum urmează: scanf("%f °/olf, &x, &xx); Notația de format %// necesită litera Z, deoarece &xx este de tip "pointer to double" Valoarea funcției returnată de funcția fscanf este egală cu numărul de elemente care au fost citite dacă au fost prezente Este egal cu zero sau cu un număr negativ dacă nu se poate citi nimic Acest fapt poate fi folosit pentru a verifica dacă s-a ajuns la sfârșitul fișierului Similar cu indicatorul stdout, există un indicator standard pentru fișierul stdin Cele două apeluri scanf( ) și fscanf(stdin, ) sunt echivalente Un singur caracter poate fi citit și scris într-un mod mai primitiv: ch = getc(fp); (intrare pe disc) putc(ch, fp); (ieșire pe disc) ' ch = getchart)', (intrare directă de la tastatură) putchar(ch); (ieșire direct pe afișaj) Aici, getcharO este exact echivalent cu getc(stdin), iar putchar(ch) este echivalent cu putc(ch, stdout) Dacă verificăm ultimul caracter citit și decidem că ar trebui folosit din nou, atunci acesta poate fi returnat înapoi la fluxul de intrare prin scriere ungetc(ch, fp)\ Funcțiile printf, fprintf, scanf, fscanf efectuează așa-numita intrare/ieșire formatate, așa cum este indicat de ultima literă f din numele lor Fișierele formatate au un liniar Anexă O SCURTĂ INTRODUCERE LA C structura și numerele din ele sunt reprezentate ca o succesiune de caractere Pentru fișierele de disc care ar trebui să fie citite și scrise numai de propriile noastre programe, va fi mai eficient să folosiți I/O neformatat Aceasta înseamnă că reprezentarea datelor internă și externă sunt identice, astfel încât numerele vor fi scrise cel mai probabil ca cuvinte binare cu lungime fixă Pentru intrare/ieșire fără format, funcțiile freadibufptr, dimensiune, n, fpY, fwritekbufptr, dimensiune, n,fp); Argumentele au următoarele tipuri și semnificații: bufptr pointer to char - pointer la bloc de memorie (numit uneori un buffer), dimensiunea int este dimensiunea în octeți a unuia element de date, n int este numărul de elemente date din buffer, fp pointer to FILE este un indicator de fișier Aceste două funcții returnează ca valoare a funcției un număr întreg egal cu numărul de elemente de date care au fost citite sau scrise Ca și în cazul funcției fscanf, poate fi folosită pentru a verifica dacă a existat o încercare de a citi date după sfârșitul fișierului, de atunci încercarea de citire va eșua și valoarea funcției fread va fi zero ALO FUNCȚII STANDARD DE MATEMATICĂ Strict vorbind, setul disponibil de funcții predefinite nu face parte din limbaj Cu toate acestea, în practică, este convenabil să aveți o listă cu astfel de funcții la îndemână Folosim adesea funcțiile matematice incluse în următoarea listă, care listează numele funcțiilor și tipurile de argumente, împreună cu o scurtă explicație despre ceea ce evaluează funcția: dublu cos(x) dublu x\ /* cosinus ♦/ dublu sin (x) dublu x\ /* sinus */ double tan (x) dublu x\ /* tangentă */ dublu log(x) dublu x\ /* logaritm natural */ A FUNCȚII MATEMATICE STANDARD dublu sqrt(x) dublu x; /* rădăcina pătrată a lui z */ etaj dublu(x) dublu x; /* cel mai apropiat număr întreg mai mic, */ /* de exemplu/Zoor( ) = */ double ceiilix) dublu x; /* cel mai apropiat număr mai mare, */ /* de exemplu celltAA) = */ int abc(i) int i; /* modul întreg */ dublu fabs(x) dublu x; /* modul virgulă mobilă */ /* punct dublu */ dublu acos(x) dublu x; /* arc cosinus */ dublu asin(x) dublu x; /* arcsinus */ dublu atan(x) dublu x; /* arc tangentă */ srandtseed) int seed-, /* inițializare generator */ /* numere aleatoare rând() */ int randi); /* generator de numere aleatorii */ long int time(p)long int p; /* timp în secunde, numărat de la ianuarie , */ /* GMT */ LITERATURĂ Ammeral, L ( ) C pentru programatori, Chichester: John Wiley & Sons Ammeral, L ( ) Grafică computerizată pentru computerul IBM, Chichester: John Wiley & Sons Ammeral, L ( ) Programe și structuri de date în C, Chichester: John Wiley & Sons Angell, I O ( ) A Practicai Introduction to Computer Graphics, Londra: Macmillan [Traducere disponibilă: Angel Y O introducere practică în grafica computerizată - M : Radio și comunicare, ] Ayres, F Jr ( ) Seria schiță a lui Schaum, Teoria și problemele geometriei proiective, New York: McGraw-Hill Escher, M C , şi colab ( ) Lumea lui M C Escher, New York: Harry N Abrams Feuer, A și N Gehani (eds) ( ) Compararea și evaluarea limbajelor de programare Ada C Pascal, Englewood Cliffs, NJ: Prentice-Hall Foley, JD și A van Dam ( ) Fundamente ale graficii interactive pe computer, citire, Mass : Addison-Wesley [Traducere disponibilă: Foley J , van Dam A Fundamentals of interactive computer graphics: In books - M : Mir, ] Forsythe, G E , M A Malcolm și C B Myer ( ) Metode computerizate pentru calcule matematice, Englewood Cliffs, NJ: Prentice-Hall Hopkins, EJ și JS Hails ( ) O introducere în geometria proiectivă plană, Oxford: The Clarendon Press Kemighan, B W și D M Ritchi ( ) Limbajul de programare C, Englewood Cliffs, NJ: Prentice-Hall [Traducere disponibilă: Kernighan B , Ritchie D Limbajul de programare C - M : Finanțe și statistică, ] Kreyszig, E ( ) Advanced Engineering Mathematics, New York: John Wiley && Sons McGregor, J și A Watt ( ) Arta graficii pe microcomputer pentru BBC Micro! Electron, Reading, Mass : Addison-Wesley " Newman, MN și R F Sproull ( ) Principiile graficii interactive pe computer, New York: McGraw-Hill [Este disponibilă o traducere a primei ediții ( ): Newman W , Sprull R Fundamentals of interactive computer graphics - M : Mir, ] Plum, T ( ) Learning to Program in C, Englewood Cliffs, NJ: Prentice-Hall index al subiectelor Clasa de memorie automată Argumente ale programului Punctul infinit , I/O Vector coloana rândul Produs încrucișat Linii verticale Vizibilitate, teste Vizualizați transformarea Vizualizați coordonatele Scara în spirală Matrice externe , Poligon convex Expresia Generarea curbei Curbe netede Ochiul Orizontul Fațete Figura chenarele Fața din spate Precizie dublă , Determinant , Diagonala Alocarea dinamică a memoriei Lungimea vectorului Fața din spate Adevărul , Dimensiunea imaginii ' Funcția pătratică constantă NULL Coordonatele lumii , ecranul Netezirea curbei Cubul , Linii verticale False , Matrice Scalare , Notație matriceală Coordonatele lumii , Observatii, con directia pct Linii ascunse Poligon neconvex I/O brut Linie nouă, caracter Normal (vector) Vector zero Zona de ieșire Dimensiunea obiectului Coordonatele uniforme , Fereastra break statement do while return sizeof , typedef Declarație de atribuire Declarații INDEX SUBIECTULUI Descriptor char float , Descriptor int static Proiecție ortografică Orificiul , în cilindrul Decuparea Parametrul argc argv epsilon variabilă Transfer , , Perspectiva Transformarea perspectivei , Pixel Piramida Arborele lui Pitagora Suprafata Virați , , direcția pozitivă Marimea Poligonul , Emisfera Sistemul de coordonate drept Transform View Perspectiva , Preprocesor , Prisma Sârmă model Program GENPLOT HIDLIN HIDLINPX , Programe, structura Geometrie proiectivă , Derivate Pseudo perspectivă Distanțe într-un triunghi Recursiune Curbe netede Limbajul de programare C Produs punctual Număr aleatoriu , Combinație de transformare Spirală Listă, liniară Spline, B-spline Corp solid , Tija Sentinel Structura Coordonate sferice , Tetraedrul Tor Punct de fugă Precizie numerică Triunghi, asociat Triunghiuri, descompunere Indicator , Expresie condiționată Robustețea programului , Dosar , stdio h Функция , Draw Endgr FABS FCLOSE , FOPEN , fprintf , FREAD FSCANF , FWRITE GETE GETCHAR INIT YWORT INITGR Main Malloc , MOVE INIT IWPORT INITGR Main Malloc , MOVE INIT IWPORT INITGR Main Malloc , MOVE INIT IWPORT INITGR Main Malloc , MOVE INIT IWPORT INITGR MAIN MALLOC , realloc scanf , srand , time ungetc Sfârșitul fișierului , I/O formatat Cilindru În sensul acelor de ceasornic, direcția de ocolire față de , Linia întreruptă , Ecranul , , Distanța până la Coordonatele ecranului Cuprins PREFAȚA Capitolul INTRODUCERE Motivație pentru nevoia de grafică programare Programare grafică în C Exercițiul Capitolul ALGORITMI D Transformare și noi coordonate împlini ani Notația matriceală Ferestre și zone de ieșire Liniile de tăiere Dimensionare și poziționare automată Utilizarea recursiunilor Curbe netede Exercițiul Capitolul UN INSTRUMENT GEOMETRIC PENTRU ALGORITMI GRAFICI D Vectorii Produs punctual Determinanți Produs încrucișat Descompunerea poligoanelor în triunghiuri Coordonate omogene' Translația și rotația în spațiul tridimensional h Exercițiul Capitolul IMAGINI DE PERSPECTIVA Introducere Vedeți transformarea Transformări de perspectivă Program de desen cub Desenarea modelelor de sârmă Direcția de observare, infinit, linii verticale Exercițiul Capitolul ELIMINAREA LINIILOR ASCUNSE Analiza performanței Date de intrare și reprezentare internă Algoritm pentru detectarea liniilor ascunse Poligoane și pixeli Programul îmbunătățit Exercițiile Capitolul EXEMPLE PRACTICE Introducere Cilindru tubular Tije spiralate Scara în spirală Thor Emisferă / Funcția a două variabile Exercițiul Aplicație O SCURTĂ INTRODUCERE LA C A Tipuri de date de bază A Unii operatori A Z Operatori și expresii A Probleme lexicale și structura programului A Matrice și pointeri A Funcții A Structuri A Alocarea dinamică a memoriei A I/O A Funcții matematice standard Literatura Index Serie GRAFICA MAȘINILOR ÎN LIMBAJ C L Ammeral (Olanda) Principii de programare în grafica computerizată Grafică de mașină pe computerele personale Grafică pe computer D interactivă Programare grafică în Turbo C (Traduceri din engleză) Patru cărți într-o formă accesibilă prezintă principalele probleme ale graficii pe computer Cărțile pot fi folosite ca ghid pentru predarea programării C cu exemple ilustrative de prezentare a rezultatelor sub formă de imagini grafice Modulele sursă ale programelor din toate cele patru cărți sunt furnizate pe dischete Include dischete de KB Pe o dischetă separată, un program de grafică D modificat pentru generarea de imagini stereo anaglife Orice formă de plată - prin transfer bancar, mandat poștal sau numerar Trimiteți comenzi pentru cărți și dischete la: Moscova K- , PO Box -Lvov Principii de programare în grafica computerizată ISBN - - -X (rusă) ISBN - - - 